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COMPUTING  INTERNAL  COCKPIT  REFLECTIONS  OF  EXTERNAL 
POINT-LIGHT  SOURCES  FOR  THE  MODEL  209  AH-1S  HELICOPTER 
FLAT-PLATE  CANOPY  DESIGN 


INTRODUCTION 

The  US  Army  Human  Engineering  Laboratory  (HEL)  has  developed  a computer  program  for 
computing  internal  cockpit  reflections  of  external  point  light  sources.  This  work  was 
accomplished  at  the  request  of  the  US  Army  Air  Mobility  Research  and  Development  Laboratory 
for  application  to  the  flat-plate  canopy  design  of  the  Model  209  AH-1S  COBRA  helicopter.  The 
flat-plate  canopy  (FPC)  design  reduces  the  solar  glint  which  occurs  during  daylight  operations  to 
a momentary  flash  at  certain  observer-aircraft-sun  angles.  Consequently,  a moving  aircraft  no 
longer  has  the  continual  solar  glint  which  occurred  with  the  earlier  compound-shaped  canopy 
designs.  The  continuous  presence  of  solar  glint  increased  the  range  of  visual  detection  by  ground 
observers. 

However,  in  certain  situations  during  night  operations  the  internal  surfaces  of  the  FPC  act 
like  mirrors  for  external  sources  of  light  as  seen  from  the  pilot’s  position.  These  virtual  images  of 
ground-level  lights  can  cause  disorientation.  The  pilot  cannot  discriminate  between  reflections 
and  light  sources  on  the  ground  and  creates  a safety  hazard  during  flight. 

The  computer  program  (Appendix  C)  was  part  of  a larger  study  effort  by  HEL  to  determine 
the  full  extent  of  the  nighttime  light  reflection  problem  and  review  various  proposals  for  solution 
(5).  The  conclusion  reached  by  HEL  is  that  the  most  realistic  solution  to  the  combined  daytime 
glint/nighttime  reflection  problem  is  a modification  of  the  FPC  side  panels  to  a simple  curved 
shape.  It  is  possible  with  little  additional  programming  to  use  the  inclosed  program  to  assist  in 
further  canopy-design  work. 


METHOD 

A three-dimensional  ray-tracing  program  was  written  to  trace  straight-line  rays  backwards 
from  the  position  (nominal)  of  the  pilot’s  eye  to  visible  points  on  the  internal  transparent 
surfaces  of  the  canopy.  Each  such  ray  is  continued  backwards  between  internal  canopy-surface 
points  until  a non-transparent  surface  is  reached.  Such  surfaces  are  assumed  to  be  diffusive 
without  specular  reflectance.  At  each  reflection  point,  the  reflectance  and  transmittance  are 
computed  along  with  the  directional  cosines  of  the  corresponding  incident  external  and  reflected 
rays.  In  this  way,  a reflected  ray  is  traced  backwards  from  the  pilot’s  eye  to  all  possible  external 
sources. 

The  method  is  applicable  to  any  canopy  design  although  application  is  limited  to  the  FPC  in 
this  report.  The  canopy  is  specified  by  a covering  net  of  flat  surfaces.  Each  surface  is  described  by 
the  coordinates  of  its  vertices  in  the  three-dimensional  space  and  the  rotational  order  in  which 
they  are  listed. 
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Tracing  rays  backwards,  the  program  computes  the  intersection  point  of  a straight-line  ray 
given  its  directional  cosines  and  origin  point,  with  a plane  containing  each  surface  segment  in 
turn.  The  reflection  point  for  the  ray  is  that  intersection  point  which  is  contained  within  the 
edges  of  the  correspt)nding  segment.  The  angle  of  incidence  between  the  surface  normal  and  ray 
at  this  point  is  computed  along  with  the  corresponding  reflectance  and  transmittance  values  and 
the  directional  cosines  of  the  transmitted  and  reflected  rays.  Working  backwards,  the  reflected 
ray  becomes  the  incident  ray  for  the  next  set  of  computations  (Appendix  A). 

The  program  includes  obstructing  surfaces  as  well  as  the  transparent  surfaces  of  the  canopy 
in  the  computations.  These  surfaces  either  obstruct  the  pilots'  forward  vision  or  block  out 
incident  rays  from  external  sources.  These  surfaces  are  (a)  the  pilot’s  display  panel  and  side 
armor,  (b)  the  copilot’s  seat,  side  armor  and  gunner’s  sight,  and  (c)  the  sides  and  floui  of  the 
cockpit.  They  are  specified  in  the  same  manner  as  are  the  canopy  segments  with  however, 
appropriate  coding  added.  The  computation  of  a reflection  point  for  a ray  on  an  obstructing 
surface  renders  the  computation  complete  since  the  backwards  traced  ray  is  considered  to  be 
absorbed. 

This  computation  process  is  repeated  for  pilot-viewing  directions  wliich  are  indexed  at  equal 
increments  over  a quarter  sector.  The  sector  is  bounded  by  vision  directly  to  the  front,  to  the 
side,  and  bottom.  In  this  way  a table  is  constructed  which  lists  at  discrete  intervals  all  possible 
internal-reflection  points  and  the  corresponding  external-light  directions.  This  approach  generates 
a large  amount  of  data  and  a computer-graphics  routine  is  included  for  output.  The  incident-ray 
entry  points  and  the  corresponding  reflection  points  are  shown  on  a perspective  drawing  of  the 
cockpit  as  seen  from  the  pilot’s  position  for  light-ray  conditions  of  interest  (sec  Appendixes  B 
and  C). 


DISCUSSION 

The  results  of  this  application  to  the  FPC  are  shown  in  Figures  1 through  4.  These  figures 
are  hard  copies  of  the  computer-graphics  output.  Figure  1 shows  side,  top,  and  front  views  of  the 
flat-plate  canopy.  Internal  obstructing  surfaces  arc  sketched  in  with  broken  lines.  The  pilot’s 
nominal  eye  position  is  shown  in  each  view.  The  aircraft  is  not  included  in  the  sketch. 

Figures  2,  3,  and  4 are  perspective  drawings  of  the  cockpit  as  seen  from  the  pilot’s  position. 
The  pilot’s  nominal  viewing  direction  is  shown  by  the  small  cross  in  the  center  of  the  upper  front 
canopy  surface.  The  figure  shows  that  the  lower  portions  of  the  front  and  forward  side  canopy 
surfaces  are  blocked  from  view  by  the  pilot’s  instrument  panel. 

The  “dots”  on  Figure  3 indicate  entry  points  of  external  rays  which  generate  visible 
reflections  on  the  right-hand  side  of  the  canopy.  The  figure  shows  that  some  entry  points  are  on 
the  lower  front  and  side  canopy  surfaces  wfiich  are  outside  of  the  pilot’s  ficld-of-view-. 

Figure  ‘I  shows  the  corresponding  reflection  points  which  arc  generated  on  the  right-hand 
side  of  the  canopy,  f he  reflection  points  are  spaced  at  2^’  x 2^  increments.  The  number  shown  at 
each  reflection  point  is  equal  to  the  negative  value  <rf  the  logarithm  (base  lO)  of  the  light 
rcllcvtanve.  i he  numlrers  are  truncated  to  the  integer  values  by  dropping  the  fractional  part.  The 
numerical  ’’/eio”  corresponds  to  those  reflectances  which  arc  greater  than  0.1  in  value.  The 
numerical  “one”  eorresponds  to  those  values  equal  to  0.1  or  less  but  greater  than  0.01 . ,\  similar 
range  exists  f<jr  the  numerical  “two.” 
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f igurc  t . Top,  side,  and  front  views  of  the  llat-pUtc  canopy  design  without  aircraft. 
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The  points  shown  on  the  figures  are  for  primary  reflectances  only.  Secondary  and  higher 
order  reflectiirns  are  possible.  Entry  rays  and  reflections  can  occur  at  surface  points  which  are 
further  back  on  the  canopy.  These  points  are  outside  of  the  perspective  view.  (See  reference  5 for 
a further  breakdown  of  entry  point  and  primary  reflection  point  distributions  by  canopy 
surfaces.) 


RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  following  recommendations  are  made  for  further  development: 

1.  Investigate  the  perspective  space  for  ray  tracing  in  place  of  the  presently  used  object 

space. 

2.  Extend  the  program  to  realistic  nighttime  lighting  situations  and  approach  scenarios. 

3.  Include  the  visual  function  of  the  pilot  as  a measure  of  the  magnitude  of  the 
reflection  problem. 

4.  Extend  the  program  to  canopy  designs  other  than  the  flat-plate  canopy. 

5.  Incorporate  similar  developments  for  the  daytime-glint  problem. 

6.  Incorporate  photometric  measurements  and  subject  testing  at  suitable  stages  in  the 
development  to  confirm  the  validity  of  the  computatory  approach. 


CONCLUSION 

A computer  program  has  been  developed  by  HEL  to  show  internal  cockpit  reflections  of 
external  point  light  sources.  The  program  has  been  applied  to  the  Model  209  AH-IS  COBRA 
helicopter  with  the  flat-plate  canopy  design.  The  results  show  that  during  nighttime  operations 
reflections  are  possible  over  much  of  the  canopy  internal  surfaces.  Which  of  tfiese  reflections 
actually  occur  depends  upon  the  particular  lighting  situation  and  flight  scenario. 
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APPENDIX  A 


INTERNAL  REFLECTIONS 


The  internal  reflections  are  computed  using  standard  vector  geometry  techniques  (1). 
Consider  a standard  cartesian  coordinate  system  (x,y,z)  with  the  y-axis  taken  along  the 
longitudinal  axis  of  the  aircraft,  the  z-axis  directed  through  the  canopy  top  and  the  x-axis 
perpendicular  to  these  two.  The  cockpit  and  canopy  may  be  described  by  a set  of  flat  planar 
surfaces.  Each  such  surface  consists  of  corners  and  straight  edges.  It  is  described  by  a sequence  of 
adjacent  corner  vertices  and  their  corresponding  coordinates,  with  the  consecutive  order 
determining  the  direction  of  the  surface  normal. 

(1)  Computation  of  surface  normal 

Consider  a flat  surface  with  three  adjacent  vertices  none  being  the  same,  P-| 
P2  Pt(x^,y:},z^).  The  surface  normal  may  be  determined  from  the 

crossproduct  and  dot  product  of  the  vectors  P^P2  and  P1P3.  These  vectors  may  be  expressed  in 
terms  of  the  unit  vectors  for  the  cartesian  coordinate  as  follows: 


P,P2  - Ax2lTx  + '^Y2i1y  + '^^21'Z 


where  AX2]  = X2  — X-|,  AY 21  = Y2  and  AZ2^  = Z2  — Let  P2  be  the 

length  of  the  vector  P]P2 

P2  = MX2,2  + AY2i2  +AZ2,2}'/2 

The  dot  product  of  the  two  vertice  vectors  is: 

P^Pz  . P^P3  = P2P3  COSO23, 

where  ©23  angle  between  the  two  vectors.  This  angle  is  measured  in  the  plane  found 

by  the  two  vectors  about  the  vertex  P^.  In  terms  of  the  coordinate  vmit  vectors, 

= AX2^  * AX3^  + AY2^  • A Y2^  + ^^21  * ^ ^31 
= ARCCOSINE  {FjP^  • 1^3} 


P,P2  . P1P3 
Therefore,  ©23 


1 1 


The  cramproduct  of  the  two  vectors  is  given  by: 


PlP2  » P1P3  = 1/2^3  "'"^^23 

v^here  tp  is  the  unit  vector  normal  to  the  surface.  In  terms  of  the  coordinate  unit  vectors, 

P^P2  • P1P3  = ( ^^21  • AZ3^  - A Z21  ‘ A Y31)  \x 
- (AX  21  ' AZ3^  - AX  3^  ' ‘^^21^  W 
+ ( AX 2^  • AY3^  - AX3^  • AY21)  I 2’ 

The  surface  normal  may  be  expressed  in  terms  of  the  coordinate  unit  vectors  as; 

Ip  = 3p  + bp  fy  + Cp  where  ap.bp,  and  Cp  are  the  appropriate 
directional  cosines.  Combining  expressions  and  using  R=l/P2P3sin  033  we  see  that 
ap  = (Ay21  ■ AX31  - Az2i  • Ay3,)  • R, 

*^n  = -(Ax2i  • AZ31  - A X3^  • A Z2j)  • R, 

Cp  = (ax2|  • Ay3i  - AX31  . Ay  21)  • R 


(2)  Specification  of  trace  ray  direction 

A straight  ray  may  be  specified  b>  an  origin  point  and  a unit  vector  directed  alorg  its 
length.  Consider  a spherical  coordinate  system  with  the  origin  at  the  pilot’s  nominal  eye  position 
and  the  y-axis  parallel  to  that  of  the  aircraft’s  cartesian  coordinate  system.  Let  a be  the  angular 
elevation  of  a ray  from  the  origin  measured  from  the  z-axis.  Lot  B be  the  azimuth  angle  of  the  ray 
neasured  from  the  x-z  plane  toward  the  y-axis.  The  direction  of  the  ray  is  specified  by  the  angles 
a , 6 . Since  the  axes  of  the  two  coordinate  systems  are  parallel,  the  directional  cosines  of  the  ray 
direction  are: 

a,;  = crjs  (B ) • sin  (a ) 

b^  = sin  ( 6 ) 

c^  = cos  (6  ) ' cos  ( a ) 
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The  origin  point  is  the  pilot’s  nominal  eye  position 


(3)  Intersection  of  ray  with  surface 

The  intersection  point  of  a line  with  a planar  surface  may  be  computed  given  the  origin 
and  directional  cosines  of  the  line  and  a vertex  and  normal  of  the  surface.  Given  the  origin  x^,  y^, 
Zj  and  the  directional  cosines  of  the  line  a5,bj,C5,  the  coordinates  of  the  intersection  point,  P,  are. 

X = X5  + a^  • R 

y = yj  + bs  ' R 
z = Z5  + C5  - R, 

where  R is  the  line  length  between  the  origin  and  the  intersection  point. 

Since  the  intersection  point,  P,  is  in  the  planar  surface,  then  the  vector  from  P to  any 
vertice  (assuming  they  are  not  the  same  point),  P^ , must  be  perpendicular ^to  the  surface  normal. 

For  this  reason,  the  dot  product  of  the  two  vectors  equals  zero,  i.e.,  P^  P-  f Therefore, 
an(x-xi)  + bn(y-y^)  + Cplz-z^)  = 0, 

and  using  the  above  expressions  for  the  intersection  point  coordinates,  the  line  length  is  given  by. 


R = {a^  • (xi-Xj)  + bp  • (yi-Ys)  + S ‘ <^T^s)  } 


*in  • ^s  + ‘’n  ' + S • S 


(4)  Determine  if  ray  intersects  surface  in  outward  direction. 

The  consecutive  order  of  vertices  assigned  for  the  cockpit  and  canopy  surfaces  is 
selected  so  that  the  surface  normals  are  directed  into  the  cockpit  interior.  The  computed 
intersection  point  will  be  outward  directed  when  the  angle  between  the  ray  vector  and  the 
surface  normal  exceeds  90°.  This  occurs  when  the  dot  product  between  the  two  vectors; 

Q = a,  . a^  + b,  • b,,  + C5  • c„. 
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is  negative,  i.e.,  Q 1 0.  Note  that  Q is  the  denominator  of  the  expression  for  R the  line  length. 
This  implies  that  for  an  outward  directed  intersection,  R^O. 

(5)  Determining  if  the  intersection  point  is  on  the  canopy  surface. 

The  intersection  point  lies  on  the  canopy  surface  when  it  is  contained  by  the  surface 
edges.  Consider  one  such  edge  defined  by  two  adjacent  vertices  Pj,Pj^-^,  and  the  ray  origin,  P^. 
These  three  points  define  a planar  constraint  surface.  A unit  normal  can  be  computed  for  this 
surface  as  above.  The  unit  normal  will  be  directed  inward  if  adjacent  points  are  selected  iri  the 
same  consecutuve  order  as  selected  for  the  canopy  surface.  The  unit  normal  for  the  constraining 
surface  is  computed  from  the  vectors  P^Pj  and  P^PjT^.  Noting  that  the  angle  between  these  two 
vectors  is  geven  by; 


j+1  = arccosine  {P^Pj  • P/j+i) 


P;  • P; 


+1 


then, 

Tc  - 

PjPj+1  sin  (0  jj+i) 


which  may  be  decomposed  into  directional  cosines  aj.,b^,c^  as  in  section  (a)  above. 

The  intersection  point  will  be  within  the  j^*^  edge  if  the  angle  between  the  line  vector 
and  the  unit  normal  of  the  corresponding  constraining  surface  is  less  than  90°.  This  is  true  if  the 
dot  product  between  the  two  vectors, 

0j  “ ^s  ■ ^cj  ■ *^s  ' *^cj  ^s  ■ ‘'cj’ 

is  greater  than  zero,  ©j>  0.  Here  j is  the  number  of  the  first  vertex  defining  the  edge.  Therefore, 
the  intersection  point  will  lie  on  a canopy  surface  with  N-vertices  if  ©j>  0 for  all  j=1 , N. 


(6)  Reflected  and  transmitted  rays 

The  canopy  material  is  assumed  to  be  dielectri^.  The  incident  ray  is  reflected  from  the 
surface  and  refracted  into  the  material.  The  unit  vector,  \^,  for  the  incident  ray  direction  may  be 
separated  into  surface  components  using  vector  notation.  The  vector  component  in  the  direction 
of  the  surface  normal  is  equal  to  -^nj^n-  while  the  component  parallel  to  the  surface  equals 
"t  tp)  Tp-  These  two  components  are  orthogonal  and  completely  define  the  vector. 

The  reflected  ray  lies  in  the  plane  defined  by  the  surface  normal  and  the  incident  ray. 
The  angle  between  the  reflected  ray  and  the  surface  normal  is  equal  in  value  but  opposite  in  sign 
to  that  between  the  incident  ray  and  the  surface  normal.  Therefore,  the  unit  vector  for  the 
reflected  ray  equals  the  surface  component  for  the  incident  ray  minus  its  normal  component,  i.e., 

I R ~"^s  " 's  ^ n^  In 
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This  may  be  expressed  in  terms  of  directional  cosines  for  the  unit  vectors  of  the  coordinate 
system, 


-2  Ran, 
-2  Rb„, 


Cfj  Cj  2 Rcp, 


where  R-a^  ' a^^b ^-bp+c^  Cp.  Note  that  the  angle  of  incidence  equals  ©Q=arccosine  (-  R). 

Consideration  of  electromagnetic  boundary  conditions  at  the  material  surface  results  in 
Ferrel’s  law.  Assuming  non-polarized  light  the  reflectance  is  equal  to  (2). 

r = (''s+'-p),  (1) 

where, 


cos  0Q  - > 

/n2  - 

sin2  9^ 

_cos  0O  + \ 

®o  - 

2 

n'^cos0Q  - 

yl^ 

- sin2 

9 

Ln  COS0  Q + 

+ sin  0^ 

The  refracted  ray  lies  also  in  the  plane  defined  by  the  incident  ray  and  surface  normal. 
The  angle  of  refraction  is  determined  by  Snell’s  law,  i.e.,  sin0Q=nsin9,  where  n is  the  material 
index  of  refraction.  The  surface  transmittance  is  determined  from  energy  considerations.  The 
energy  flux  in  the  ray  beam  is  equal  to  the  product  of  the  intensity  (i.e.,  square  of  electric 
amplitude  times  dielectric  constant)  and  the  beam  crosssectional  area.  The  conservation  of  energy 
specifies  that  the  sum  of  the  energy  fluxes  in  the  reflected  and  refracted  beams  equals  that  in  the 
incident  beam.  Assuming  distant  sources  of  light  so  that  the  wavefronts  are  closely  planar,  the 
surface  reflectance  and  transmittance  is  related  by; 


r+t 


cos  9 = 


cos 


% 


1 


(2) 


The  refracted  ray  is  multiply  reflected  between  the  top  and  bottom  surfaces  of  the 
canopy  material.  This  internally  transmitted  energy  is  partially  absorbed  as  it  transverses  the 
medium  between  the  surfaces.  At  each  surface  the  energy  remaining  is  partially  refracted  from 
the  medium  and  the  rest  reflected  toward  the  other  side.  The  same  reflection  and  refraction 
equations  apply  for  this  internal  reflections  as  above  except  that  the  angles  0^  and  0 are 
interchanged. 

The  total  reflectance  and  transmittance  can  be  determined  by  adding  together  the 
fractions  of  energy  that  leave  the  top  and  bottom  surfaces.  Assume  that  the  top  and  bottom 
surfaces  are  parallel  and  that  the  material  thickness  is  infinitesimal.  All  internally  refracted  rays 
leave  at  the  same  surface  point  as  does  the  initially  transmitted  ray,  and  this  point  is  the  same  as 
that  at  which  the  incident  ray  reaches  the  surface. 

Let  r^  and  t^  be  the  surface  reflectance  and  transmittance  for  the  reflection  and 
refraction  of  an  externally  incident  ray  (see  equations  1 and  2).  Similarly,  let  rj  and  tj  be  the 
internal  surface  reflectance  and  transmittance  (equations  1 and  2 with  0 substituted  for  0,j). 
Letting  T j be  the  internal  transmittance  of  the  material,  the  total  reflectance  is  given  by  (3). 


= rg  + ty  ■ T|  (rj  ■ Tj  tj  +rp  Tj^  rj  • tj  ■ +....) 


and  the  total  transmittance  by 


T = t^Tj(tj-  +R.2.t2.  tj-  + 


In  closed  form,  these  equations  are 


(1-r?T? 


T='otjTj 


(1-r?T2  ) 


The  total  transmitted  ray  will  exit  at  the  intersection  point  and  have  directional  cosines  equal  to 
those  of  the  incident  ray,  i.c.,  aj=a5,bj=b^,  and  cj=c^. 


(7)  Repeated  computations  and  interpretation 


The  above  computation  process  can  now  be  repeated  to  determine  a canopy 
intersection  point  for  the  reflected  ray.  Put  the  incident  ray  equal  to  the  just  computed  reflected 
ray,  i.e.,  ‘is=a|^,bj=b[^,C5=C|^,  and  the  ray  origin  to  the  intersection  point,  i.e., 

X5=X|,Y^=Y|,Z^=Z|.  Continuing  the  computation  process  will  determine  the  next  intersection 
point,  the  corresponding  reflected  and  transmitted  rays,  and  the  reflectance  and  transmitted 
values. 

Assume  that  N such  successive  reflections  have  been  determined  starting  with  a ray 
from  the  pilot’s  nominal  eye  position.  One  interpretation  of  the  results  is  that  the  transmitted  ray 
of  the  N—  computation  stage  is  the  negative  of  an  incident  ray  from  an  external  light  source. 
This  ray  is  transmitted  through  the  canopy  at  the  intersection  point  with  transmission  T(N). 
The  ray  is  then  reflected  in  the  direction  opposite  to  that  computed  until  it  reaches  the  pilot’s 
eye  position.  The  final  reflectance  value  is  equal  to: 

R=T(N)-  R(i), 


where  R (i)  is  the  reflectance  computed  at  the  i—  intersection  point.  Note  that  for  N > 1, 
R(0)=1. 


APPENDIX  B 

PERSPECTIVE  DRAWINGS 


Perspective  drawings  of  the  cockpit  as  seen  from  the  pilots’  nominal  eye  position  have  been 
prepared  using  a computer-graphics  terminal.  The  pilots’  viewing  direction  is  assumed  to  be 
fixated  along  a line  which  is  directed  five  degrees  above  the  aircraft’s  longitudinal  axis  (negative 
y-axis).  The  fixation  linelies  in  the  y-z  plane  of  the  aircraft’s  coordinate  system. 

A perspective  drawing  shows  a three-dimensioal  view  of  objects  on  a two-dimensional  scene. 
The  drawing  should  be  held  by  the  viewer  at  a specified  distance  (say  a)  from  his  eye  position  and 
in  a plane  normal  to  this  fixation  direction.  The  two-dimensional  scene  is  normally  bordered  on 
four  sides  (each  of  length,  say  2b).  Not  shown  on  the  scene  are  objects  which  would  appear 
outside  of  this  border  or  behind  the  viewer.  The  scene  may  be  prepared  on  a graphics  terminal  of 
dimension,  say  2c.  The  drawings  shown  in  this  report  fill  a 10-inch  by  10-inch  scene  and  have  a 
10-inch  viewing  distance.  The  graphics  terminal  used  in  preparation  has  a span  of  1024  rasters 
(a=10  inches,  b=5  inches,  and  c=512). 

The  viewer’s  eye  position  may  be  considered  as  forming  the  apex  for  a four-sided  viewing 
pyramid.  The  sides  of  the  pyramid  are  determined  by  the  apex  and  the  borders  of  the  scene.  The 
objects  within  the  viewing  pyramid  are  projected  as  images  onto  the  scene  surface.  A straight  line 
connecting  the  apex  with  a vertex  of  an  object  will  pass  through  the  corresponding  vertex  of  the 
object’s  image.  Consequently,  the  image  is  described  by  its  vertices  in  the  scene  and  the  edges 
connecting  them. 

Each  object  is  described  by  its  surfaces  and  the  corresponding  vertices  and  edges.  Only  the 
images  of  surfaces  which  are  facing  the  viewer  should  be  drawn,  since  the  back  surfaces  are 
hidden  from  view.  In  some  cases,  a surface  will  be  partially  hidden  from  view  by  intermediate 
objects  and  only  an  image  of  the  visible  portion  should  be  drawn.  Still  other  objects  may  be 
completely  hidden  from  view  by  intermediate  objects.  Other  objects  will  be  outside  of  the 
viewing  pyramid  and  completely  hidden  from  view.  Finally,  some  remaining  objects  will  be 
partially  inside  the  viewing  pyramid  and  partially  outside. 

The  construction  of  perspective  drawings  follows  a three-stage  decision  process:  First,  find 
those  surfaces  which  face  toward  the  viewer.  Next,  find  which  of  these  surfaces  arc  partially  or 
completely  within  the  viewing  pyramid.  Finally,  determine  what  portions,  if  any,  of  the  surface 
edges  are  not  hidden  by  the  sides  of  the  pyramid  or  intermediate  objects. 

The  above  decision  process  is  implemented  by  a seven  step  analysis  sequence  (4). 

1 . Compare  each  surface  normal  to  the  viewer’s  direction  and  discard  surfaces  which  do 
not  face  the  viewer. 

2.  Convert  the  coordinates  of  the  surface  vertices  from  the  aircraft  object  space  to  the 
viewer’s  space. 

3.  Trim  off  surface  edges  which  extend  behind  the  viewer. 
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4.  Convert  the  vertex  coordinates  from  the  viewer’s  space  to  a scene  display  space  in  j 

which  the  viewing  pyramid  has  become  a view  box.  j 

5.  Discard  surfaces  which  lie  completely  outside  the  view  box.  ] 

6.  Determine  what  portion  of  a surface  edge  lies  within  the  viewing  box.  | 

7.  Determine  the  apparent  crossing  points  of  the  surface  edge  with  the  edges  of  every  j 

other  surface.  Determine  whether  the  points  are  on  constraining  edges  of  another  surface,  I 

whether  they  are  in  front  of  the  test  edge,  and  finally,  what  portion  is  blocked  by  the  , 

constraining  surface.  i 

This  sequence  determines  the  end  points  of  the  visible  portion  of  each  surface  edge  and  the  l 

resulting  image  lines  are  drawn  on  the  perspective  drawing. 

(1)  Discarding  surfaces  not  facing  the  viewer. 

The  consecutive  ordering  of  surface  vertices  was  chosen  so  that  surface  normals  are  are 
directed  into  the  cockpit  volume.  Consider  a line  from  the  pilot’s  nominal  eye  position  (xQ,yQZQ) 
to  a vertex  on  a surface  (x^  ,y  ^ ).  The  directional  cosines  (a5,b^,Cj)  of  a unit  vector  directed 

along  this  line  are  computed  from  the  differences  in  coordinate  values  and  the  line  length,  R^, 
i.e., 

a^  — (x^-Xq)  /Rj, 

= (Vl  "Vo) 

Cj  - (^1  - Zy)  /Rj 

I 

Consider  now  the  dot  product  between  the  unit  vector  along  this  line  and  the  unit  normal  at  the 
surface.  The  unit  normal  is  described  by  its  directional  cosines  (af^,b,.,,c^)  and  the  dot  product  is 
given  by  | 

e = as  ■ an  + bj  • bn  + C,  • c^ 

This  value  is  equal  to  the  cosine  of  the  angle  betwen  the  two  vectors.  The  surface  normal  is 
directed  back  toward  the  viewer  when  this  angle  is  greater  than  90*^.  In  this  case  the  dot  product 

value  is  less  than  zero,  i.e.,  0 < 0.  j 

(2)  Converting  the  coordinates  of  surface  vertices  to  the  viewer’s  space.  | 

The  coordinates  of  the  surface  vertices  are  transformed  from  the  aircraft’s  coordinate  | 

system  (x,y,z)  to  the  viewer’s  coordinate  system  (x  ,y  ,z‘)  in  two  steps.  The  origin  of  the 
coordinate  system  is  translated  to  the  pilot’s  nominal  eye  position.  The  coordinate  axes  are  then 
rotated  about  the  origin  so  that  the  transformed  y -axis  is  aligned  along  the  pilot’s  nominal 
fixation  direction.  The  x -axis  remains  in  the  x-y  plane  of  the  aircraft  space. 

I 

i 
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The  aircraft’s  coordinate  system  is  a left-handed  cartesian  coordinate  system  with  the 
y-axis  directed  backward  along  the  aircraft’s  longitudinal  axis.  Consider  now  a translation  of  the 
coordinate  origin  to  the  pilot’s  nominal  eye  position,  and  let  the  y-axis  be  directed  in  the 
opposite  direction  to  form  a right-handed  cartesian  coordinate  system.  The  pilot’s  fixation 
direction  may  be  specified  in  this  coordinate  system  by  the  spherical  angles  a and  (see 
Appendix  2A)).  The  angular  elevation,  a,  of  a ray  from  the  origin  is  measured  from  the  z-axis. 
The  azimuth  angle  3 is  measured  from  the  x-z  plane  in  the  x-y  plane. 

The  coordinate  system  may  be  aligned  with  the  fixation  direction  by:  (1)  a clockwise 
rotation  about  the  z-axis  equal  to  the  complement  of  the  angle,  3 , and  (2)  a counterclockwise 
rotation  about  the  x-axis  equal  to  the  complement  of  the  angle,  a-  The  rotations  are  performed 
about  the  origin  in  the  appropriate  directions  when  viewed  from  the  positive  sides  of  the  axes. 
The  above  translational  and  rotational  transformations  may  be  described  by  concatenated 
matrices. 

Let  (x,y,z)  be  a pointin  the  aircraft’s  space  and  (x^y^.z^)  the  equivalent  point  in  the 
pilots’  viewing  space,  then  in  matrix  notation  — 

(x^y^z^1)  = (x,y,z,1)  • Mj  • • Mx 


The  translational  matrix  is  given  by; 


My  = 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

-Yo 

-^o 

1 

and  the  two  rotational  matrices  by: 

Mz  = 

sin  3 

- cos3 

0 

0 

cosB 

sin  3 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

and 

Mx  = 

1 

0 

0 

0 

0 

sin  a 

cosa 

0 

0 

-cosx 

sin  a 

0 

0 

0 

0 

1 

- 

- 

, 
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The  order  of  application  of  the  transformations  are  to  be  preserved  when  the  matrices 
are  multiplied  together.  The  resulting  equations  are; 

X’  = X • sing  - (y-yjcosg 

yl  = - X • cosgsina-  (y-yJsingsinDi  - (z-z^jcosa 
= -xcosgcosa-  (y-yQ)singcosa  + (z-z^lsina 

In  the  case  considered  in  this  report,  g =90°  and  a =85°.  The  above  equations  reduce  to: 

X^  = X 

= -(y-yo)cos(5°)  - (z-zjsin(5°) 

Z’  = -(y-Yo)sin(5°)  + (z-Zo)cos(5°) 


(3)  Trimming  surface  edges  lying  behind  the  viewer  — 

Our  objective  is  to  place  a two-dimensional  scene  on  a plane  which  is  normal  to  the 
fixation  direction  (i.e.,  the  viewer’s  y-axis)  and  at  a distance  a from  the  viewer's  eye.  The  scene 
is  to  be  identical  in  appearance  to  the  three-dimensional  view  contained  within  the  scene  borders. 
It  will  be  to  our  advantage  to  trim  surface-edges  which  extend  behind  the  viewer,  for  reasons 
given  later. 

Consider  a trimming  plane  normal  to  the  y-axis  and  placed  a minute  distance,  y^,  in 
front  of  the  eye  position.  Surfaces  behind  the  trimming  plane  are  to  be  discarded,  while  surfaces 
penetrating  the  plane  are  to  be  trimmed  to  reduced  forms.  In  a sense,  we  are  considering  a 
truncated  viewing  pyramid,  since  surfaces  between  the  trimming  plane  and  the  eye  position  will 
be  discarded. 

Consider  the  starting  and  ending  vertices  of  each  edge  on  a surface.  If  the  value  of  the 
y-coordinate  of  a vertex  point  is  less  than  y^,  then  the  vertex  must  be  behind  the  trimming  plane. 
Let  the  edge  be  defined  by  the  vertex,  Pj(Xj,Yj,Zj),  for  a starting  point  and  the  j+1^^ 
vertex,  Pj+^(Xj+pYj+^  ,Zj+^ ),  for  an  ending  point.  The  following  cases  are  possible;  (1)  if  yj  and 
Yj+i  are  both  larger  than  y^,  then  the  j—  edge  is  inside  the  trimming  edge,  (2)  if  Yj  > y^,  but  yj+^ 
< yj,  then  the  edge  crosses  out  of  the  trimming  plane,  (3)  if  yj  and  yj+^  < y^,  then  the  edge  lies 
outside  the  plane,  and  (4)  if  yj  < y^,  but  yj+^  Z y^.,  then  the  edge  crosses  inside  the  plane. 

We  may  construct  an  equivalent  but  reduced  surface  by  defining  a set  of  vertices  (Pgi^} 
as  follows: 

(1)  If  yj  and  yj+^  are  both  larger  than  y^,  then  set  the  coordinates  of  P^|^  equal  to 
those  of  Pj  i.e.,  Xek=Xj,ye|,=yj,  z^|,=Zj. 


(2)  If  ^ Z yj,  but  yj^.]  < y^,  then  the  coordinates  of  P^j^  equal  to  those  of  Pj  as 
[)sc  of  PpL-  + i equal  to  those  of  P|;+i  • Here  Pi  is  the  intersection  point  of  the  cdj 


but  set  those  of  Peg+t  equal  t( 
the  trimming  plane  (see  below). 


as  above, 
edge  with 


(3)  If  Yj  and  y■^+-^  < y,,  then  discard  the  edge. 

(4)  If  Yj  < Yf  but  yj+)  > Yp  then  set  the  coordinates  of  P^|.^  equal  to  those  of  P| 
defined  above. 

Note  that  if  all  vertices  arc  outside  of  the  trimming  plane,  this  construction  will  discard  the 
corresponding  surface. 

The  coordinates  of  the  intersection  point  that  an  edge  makes  with  the  trimming  plane 
may  be  computed  using  the  technique  of  Appendix  3A.  Note  that  the  directional  cosines  of  the 
normal  to  the  trimming  plane  are  (0,1,0)  and  that  a convient  vertex  point  for  this  plane  is 
(0,Yj,0).  Consider  the  edge  with  starting  point,  Pj  and  ending  point,  Pj+i  - The  directional 

cosines  for  a unit  vector  along  the  edge  are  given  by; 

a^  — (xj^.^  — Xj)/Rg, 

bj  = (Yj+i-yj)/Re> 

^s  ~ (^j+1  ~ 

where  the  edge  length  is  given  by: 

Re  = ^ (yj+i-y,)^  + 

The  distance  along  the  edge  from  the  starting  point  to  the  intersection  point  is  given  by  (see 
Appendix  3A). 

R|  = (yt-yj)/bs- 

The  coordinates  of  the  intersection  point  are  given  by; 

X|  = Xj  + a^R, 

Y,  = Yj  + b,R, 
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(4)  Converting  the  coordinates  of  surface  vertices  to  a scene  display  space. 

The  scene  coordinates  (Xq,  Yq)  for  an  image  point  of  a surface  vertex,  P{x,y,z)  in  the 
viewer’s  space,  are  defined  by  the  intersection  of  a straight  line  with  the  scene  plane.  The  line 
runs  from  the  viewer’s  eye  position  to  the  vertex  point.  This  line  along  with  the  y-axis  of  the 
viewer’s  space  forms  two  congruent  triangles,  one  for  the  image  point  and  the  other  for  the 
vertex  point.  The  ratio  of  the  eye-scene  distance,  a,  to  the  image-point  distance  normal  to  the 
y-axis,  is  equal  to  the  ratio  of  the  y-distance  for  the  vertex  point  to  the  corresponding  normal 
distance,  i.e. 

Xq  = 2L  ■ 3, 

Y 

Zq  = Z_  ■ a. 

Y 

Consider  the  addition  of  a third  coordinate  to  add  depth  thereby  forming  a display 

space, 


Yd  = -1/Y. 

This  transformation  converts  the  viewer’s  eye  position  to  minus  infinity.  All  lines  from  the  eye 
position  to  surface  vertex  points  are  parallel  to  the  Yp-axis.  The  (truncated)  viewing  pyramid  is 
converted  into  an  infinitely  deep  viewing  box.  The  four  sides  of  the  viewbox  parallel  the  Yp-axis. 
Assuming  that  the  scene  borders  parallel  the  scene  coordinate  axes  (Xp,Yp),  then  the  viewbox 
sides  lie  in  the  planes  Xp=+b  and  Zp=+b. 

Edges  which  pass  through  the  trimming  plane  in  the  viewer’s  space  are  difficult  to  trace 
in  the  display  space.  The  portion  of  the  edge  in  front  of  the  plane  is  transformed  into  a line 
running  from  the  endpoint  image  to  positive  infinity.  On  the  other  hand,  the  portion  behind  the 
plane  is  transformed  into  a line  running  from  negative  infinity  to  the  other  endpoint  image.  This 
is  not  readily  apparent  when  the  two  endpoints  are  considered  alone.  Note  also  that  placing  the 
trimming  plane  at  a minute  distance,  Y^,  in  front  of  the  viewer  precludes  divisionby  zero  when 
transforming  trimmed  edges. 

(5)  Discarding  surfaces  completely  outside  of  the  display  space  viewbox. 

A surface  vertex  will  not  be  seen  if  the  image  coordinates  exceed  the  scene  bounds  of 
length  2b.  Consider  the  coordinate  conversions; 

^ ^D/b>  ^ ^D/b- 


If  either  the  absolute  value 
image  should  not  be  in  the  scene. 


ofXp  or  Zp, 


or  both  are  larger  than  unity,  then  the  vertex 
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(6)  The  clipping  of  surface  edges  protruding  beyond  the  viewing  box  — 

The  surface  edge  is  tested  against  the  viewbox  in  display  space.  A starting  point 

directional  cosines  (ay,by,Cy)  and  edge  length,  R^,  are  computed  from  the  end  vertices 
(Pj,Pj+l ).  Let  Rq  be  the  distance  along  the  edge  line  from  the  starting  point  to  the  initial  point  of 
the  edge  portion  within  the  viewbox  and  Rp  the  distance  to  the  final  point.  Initially  set  RO=0 
and  Rp=R^  when  entering  the  testing  procedure.  Note  that  if  Rq  is  larger  than  Rp,  then  no 
portion  of  the  edge  is  within  the  viewbox. 

The  unit  vector  normal  to  the  surface  of  a viewbox  side  is  assumed  to  be  directed  into 
the  viewbox  interior.  Letting  the  directional  cosines  of  the  surface  normal  be  denoted  by 
a^,bp,Cp,  a convenient  vertex  (x^,yp,,£p)  within  the  surface  side  is  determined  by  Xj^=-b-ap,, 
y„-b-b„,andZn=-b-Cn. 

Consider  the  dot  product  of  the  unit  vector  taken  along  the  edge  line  and  the  surface 
normal  of  the  viewbox  side,  i.e., 

Q - ^v  ■ 3n  + by  ■ bn  + Cy  • Cn 

If  the  vectors  are  orthogonal,  then  the  angle  between  them  equals  90°  and  Q equals  zero.  In  this 
case,  the  edge  line  parallels  the  side  surface,  and  the  starting  point  (xy,yy,Zy)  must  be  tested  to 
determine  whether  or  not  it  is  within  the  viewbox  (see  section  5,  above). 

If  the  dot  product,  Q,  is  not  equal  to  zero  then  the  edge  line  intersects  the  side  surface. 
The  distance  along  the  edge  line  from  the  starting  point  to  the  intersection  point  is  given  by  (see 
Appendix  3A). 

R|  = ( (X|-Xy)  • a„  + (y„-yv)  • bn  + (z^-Zy)  • c^)IQ. 


If  the  edge  line  crosses  into  the  side  surface  from  outside  the  viewbox,  then  the  angle 
between  the  two  vectors  is  less  than  90°  and  the  value  of  Q is  greater  than  zero.  In  this  case  given 
that  the  distance  R|  is  greater  than  R^,  then  the  initial  point  of  the  contained  portion  must  be 
moved  forward  to  the  intersection  point,  i.e.,  Rq=R|. 

If,  however,  the  edge  line  is  directed  outward  from  inside  the  viewbox,  then  the  angle 
between  the  two  vectors  is  greater  than  90°  and  Q is  less  than  zero.  In  this  case  given  that  the 
distance  i<|  is  less  than  Rp  then  the  end  point  of  the  contained  portion  must  be  moved  back  to 
the  intersection  point,  i.e.,  Rp=Rp 
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(7)  Removal  of  hidden  lines  on  surface  edges  and  display  of  visible  segments. 

The  clipping  routine  determines  which  segment  of  a surface  edge  is  contained  within 
the  viewbox.  The  test  edge  is  specified  by;  (1)  a starting  point  (x^.y^/^),  (2)  directional  cosines 
(a^.by.c^)  and  (3)  a length,  R^.  The  contained  segment  is  determined  by  its  initial  and  final 
endpoints  on  the  edge.  The  endpoints  can  be  described  by  their  distances  along  the  edge  line 
from  the  starting  point.  Let  R^  denote  the  distance  of  the  initial  point  and  Rp  that  of  the  final 
point.  Note  that  if  R^  equals  ^ero  and  Rp  equals  R^.  than  the  complete  edge  is  contained  within 
the  viewbox.  If  R^  is  greater  than  zero  or  Rp  is  less  than  R^  or  both,  then  only  a portion  of  the 
edge  is  contained.  However,  if  Ry  exceeds  Rp  in  value,  than  no  segment  can  be  contained. 

The  contained  segment  of  the  edge  may  be  cither  partially  or  completely  hidden  from 
view  by  intermediate  surfaces.  Assume  all  surfaces  to  be  convex.  That  is,  any  point  defined  as  the 
.linear  sum  of  edge  points  each  times  a coefficient  the  sum  of  which  equals  unity,  is  on  the  surface 
contained  by  the  edges.  In  this  case,  a line  running  along  the  test  edge  will  appear  to  the  viewer  to 
intersect  a constraining  surface  at  two  points.  Each  point  is  the  apparent  crossover  of  the  test  line 
with  an  edge  of  the  constraining  syrfacc.  For  this  reason,  an  intermediate  surface  can  mask  no 
more  than  a single  line  segment  from  view. 

Consider  the  apparent  crossover  point  between  the  test  edge  and  a constraining  edge.  A 
starting  point  (Xg,yg,Ze),  directional  cosines  (ap,b^,,Cg)  and  line  length  R^,  may  be  computed  for 
the  constraining  edge  given  the  coordinates  of  the  starting  and  ending  vertices  (Pj,  Pj+i).  The 
image  lines  of  the  two  edges  in  the  display  scene  will  not  appear  parallel  if  the  y-component  of 
the  crossproduct  is  unequal  to  zero,  i.e.,  the  value  of  Q^a^-  c^-a^-c^,,  does  not  equal  zero.  In  this 
case  there  is  an  apparent  crossover  point.  It  is  determined  in  the  display  space  by  the 
intersections  of  the  two  edges  with  a straight  line  parallel  to  the  y-axis  and  from  the  viewer’s 
position.  The  intersection  points  have  the  same  display  coordinates  (x|,y|),  but  different  dept 
values  (yie.yiv)-  *^lv  denote  the  distance  from  the  starting  point  of  the  test  edge  to  the 
intersection  point  with  the  viewing  line.  Let  R|g  denote  the  equivalent  distance  for  the 
constraining  edge.  The  display  coordinates  for  the  intersection  points  are  given  for  the 
constraining  edge  by  — 

X|  - Xe  ae  ■ R|e> 

^1  “ ^e  ^ '"e  ‘ ^ le’ 
and  for  the  test  edge  by: 

X|  = Xv  + a^  . R|y, 

Z|  = Yv  + S ■ 


Solving  the  above  equations  for  the  intersection  distance  on  the  test  edge  results  in: 


^Iv  ^^e  ' ^^e  ^v^  ^e  * ^^e  ^v^  )/Q- 

Similarly,  solving  for  the  intersection  distance  on  the  constraining  edge  results  in: 
R|e  = (^v'^e  'V  • '^lv)/‘^e- 
or  equivalently, 

R|e  = (Xy-Xg  + ay  . RiJ/a^. 


The  constraining  edge  will  be  part  of  a masking  surface  if  two  conditions  exist.  First, 
the  intersection  point  must  be  within  the  edge  endpoints.  That  is,  the  distance  R|j.  must  be 
greater  than  zero,  i.e.,  R|g  > 0,  and  less  than  the  edge  length,  i.e.,  Rj^  < Rg.  Furthermore,  the 
constraining  edge  must  be  between  the  viewer  and  the  test  edge.  The  distance  along  the  viewing 
direction  (i.e.,  the  y-axis)  to  the  intersecting  point  on  the  test  edge  is  given  by: 

Vlv  = Vv  + ^ • "^Iv- 

and  the  distance  to  the  intersection  point  on  the  constraining  edge  is  given  by 
Vie  “ ^e  ^ *^e  *^le' 

The  constraining  edge  will  appear  in  front  of  the  test  edge  if  the  distance  V|g  is  larger  than  y |y, 

'•e-yie^yiv 


A masking  surface  will  intersect  the  test  edge  line  at  two  points  Let  the  distance  trrim 
the  test  edge  starting  point  to  the  closest  intersection  point  be  denoted  b\  R^,  while  that  to  ttie 
farthest  intersection  point  be  R-^.  Note  that  R2  is  larger  than  Rj,  i.e.  R2>  R].  Four  masking 
cases  are  possible: 

(1)  The  surface  may  mask  a segment  of  the  test  edge  line  which  is  outside  of  the  edge 
limits.  In  this  case  Rj  is  larger  than  the  edge  length,  R^  or  R2  is  less  than  zero.  The  constraining 
surface  does  not  mask  the  test  edge  and  may  be  ignored. 

(2)  The  surface  may  mask  the  entire  edge.  In  this  case  R]  is  less  than  zero  and  R2  is 
greater  than  Ry.  No  segment  of  the  edge  is  visible. 
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(3)  The  surface  may  mask  a central  portion  of  the  test  edge.  In  this  case  is  greater 
than  ^ero  and  R2  is  smaller  than  R^.  The  masked  segment  extends  from  R^  to  R2. 

(4)  The  surface  masks  one  end  point  of  the  edge  but  not  the  other.  Here  the  above 
three  conditions  do  not  apply,  and  (a)  R^  is  greater  than  zero  or  (b)  R2  is  less  than  R^.  In  the 
former  case  the  masked  segment  extends  from  the  starting  point  to  R-|.  In  the  latter  case  it 
extends  from  R2  to  the  end  point. 

The  test  edge  may  be  masked  by  more  than  one  constraining  surface.  Consider  a file  of 
N visible  segments  tor  the  test  edge,  { RO(i  ),RF(i  ),  i = l ,N  } . Here  RO(i  ) is  the  distance  from  the 
edge  starting  point  to  the  initial  point  of  the  i^  segment  and  RF(i  ) is  the  corresponding 
distance  for  the  end  point.  Initially,  this  file  may  be  generated  using  the  end  point  distances  for 
the  edge  segment  contained  within  the  viewbox,  i.e.,  N=1 , RO(1  )=Rq  and  RF(1  )=Rp  (see  section 
6 above).  The  test  edge  may  then  be  compared  against  each  of  the  other  surfaces  in  turn  and  the 
change  in  the  display  file  computed  as  each  masking  surface  is  located. 

Suppose  that  a masking  surface  is  in  one  of  the  above  cases.  The  corresponding  change 
to  the  display  file  is  as  follows; 

(1)  The  constraining  surface  does  not  mask  the  test  edge.  No  change  in  visible  segments 

occurs. 

(2)  The  constraining  surface  masks  the  entire  edge.  The  number  of  visible  segments 
wihtin  the  file  is  set  equal  to  zero. 

(3)  The  constraining  surface  masks  a central  portion  of  the  test  edge  which  extends 

from  R^  to  R2.  Considering  each  of  the  visible  segments  in  turn  we  note  that  for  the  i — 

segment: 

(a)  If  RF(  i ) is  greater  than  R^  but  less  than  R2  the  segment  protrudes  into  the  masked 
region  and  RF(i  ) must  be  set  equal  to  R^. 

(b)  If  RO(  i ) is  greater  than  R^  and  RF(  i ) is  less  than  R-j,  the  segment  is  masked  and  it 
must  be  removed  from  the  file  and  the  file  number  reduced  by  one. 

(c)  If  RF(i)  is  greater  than  R2  and  RO(i  ) is  less  than  R^,  the  segment  straddles  the 

masked  section.  The  segment  is  separated  into  two  segments  by  setting  RO(k)  = R2,RF(k)=RF(i  ) 
and  RF(i  )=R^.  The  file  number  is  increased  by  one  to  include  the  segment,  k. 

(d)  If  RO(i ) is  greater  than  R^  but  less  than  R2  and  RF(i  ) is  greater  than  R2,  the 

segment  protrudes  out  of  the  masked  region.  In  this  case  RO(i  )=R^. 
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(4)  The  constraining  surface  masks  one  end  of  the  test  edge.  If  the  starting  end  is 
masked  let  Ryo~*^2  terminal  end  is  masked  let  Ry^^O  and  RypR^.  Then 

considering  each  of  the  visible  segments  in  turn,  we  note  that  for  the  i^^  segment: 

(a)  If  RF(  i ) is  less  than  Ry^  the  segment  is  masked;  it  must  be  removed  from  the  file, 
and  the  file  number  reduced  by  one. 

(b)  If  RO(i)  is  less  than  the  segment  protrudes  out  of  the  masked  region  and 
RO(i)=RyQ. 

(c)  If  RF(i)  is  greater  than  Ryf,  the  segment  protrudes  into  the  masked  region  and 
RF(i)=Ryf. 

(d)  If  RO(i ) is  greater  than  Ry^  the  segment  is  masked;  it  must  be  removed  from  the 
file  number  reduced  by  one. 

The  result  is  a file  of  N visible  segments  for  each  surface  edge.  The  file  contains  an 
initial  and  final  distance  for  each  visible  segment.  The  scene  coordinates  for  the  image  points  of 
each  of  the  segment  vertices  are  for  the  initial  vertex, 

X|  = Xy  + 3y  • RO  ( i ), 

= z + c • RO  ( i ), 

V 

and  for  the  final  vertex, 

xp  = Xy  + ay  • RF  ( i ) 

zp  = Zy  + by  ■ RF  (i  ). 


The  scene  coordinates  are  converted  to  the  display  screen  coordinates  (Xj,/.^)  using  x^=(xq+1)-  C 
and  Z5=(zq+1.  Here  2c  is  the  screen  dimension  (C=512).  A display  line  element  is  then  drawn 
between  the  starting  and  ending  vertices  for  each  of  the  N visible  segments. 
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APPENDIX  C 


COMPUTER  PROGRAM 


The  computer  program  is  programmed  for  a disk-based  batch  operating  system  in  Fortran 
IV  language  (6).  A listing  of  the  subroutines  follows.  The  program  is  attached. 

1.  CONTE  - controls  canopy  data  read  in,  computation  of  reflection  points,  ray 
indexing,  file  writing,  and  data  output. 

2.  READV  - reads  in  from  cards  the  number  of  vertices,  the  number  of  canopy 
surfaces,  the  coordinates  of  each  vertex,  the  number  of  vertices  for  each  surface,  and  the 
assignment  of  vertices  to  each  surface  in  rotational  order.  Reads  in  also  coordinates  of  pilot’s 
nominal  eye  position.  Called  by  CONTE. 

3.  NORML  - uses  coordinates  of  vertices  and  assignment  of  vertices  to  canopy  surfaces 
to  determine  coordinates  of  surface  vertices.  Uses  these  coordinates  and  order  of  vertice 
assignment  to  compute  surface  normals.  Writes  computed  data  into  disk  file.  Called  by  CONTE. 

4.  READF  - reads  in  previously  stored  canopy  data  from  desk  file  to  common  area. 
Called  by  CONTE. 

5.  TEETY  - prints  out  canopy  data  from  common  area.  Called  by  CONTE. 

6.  INTEC  - computes  intersection  point  of  straight  line  with  planar  surface  given 
origin  and  directional  cosines  of  line  and  coordinates  of  surface  vertex  and  directional  cosines  of 
surface  normal.  Determines  if  intersection  point  is  contained  by  edges  defined  by  surface  vertices. 
Called  by  CONTE. 

7.  COMP  - computes  incident  angle  between  ray  and  surface  normal  at  intersection 
point,  and  the  corresponding  reflectance  and  transmittance.  Called  by  CONTE  as  a function  of 
INTEC  return. 

ASIN 

ACOS  - computes  arcsine  and  arccosine  of  function  value.  Called  by  COMP  and  other 

routines. 

9.  DRWCN  - generates  side,  top  and  front  views  of  canopy  as  a check  on  data  input. 
Initiates  display,  sets  up  display  entities  in  display  file  using  canopy  data  in  common  area,  and 
starts  display.  Uses  Fortran  graphics  commands.  Called  by  CONTE. 
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10.  EMAKK 


ETIC  - Loads  display  tile  with  line  elements  for  display  entity  showing  nominal  eye 
position  in  canopy  views.  Called  by  DRWCN.  ETIC  called  also  by  PERP. 

1 I.  LINS 

LINE 

LINT  - loads  display  file  with  line  elements  for  canopy  surface  entities.  Called  by 

DRWCN. 


12.  PERP  - generates  perspective  drawing  of  cockpit  as  seen  from  pilot’s  nominal  eye 
position  and  shows  entry  and  reflection  points.  Initiates  display,  loads  display  file  using  canopy 
data  in  common,  reads  disk  file  for  computed  entry  and  reflection  data,  selects  data  for  showing 
according  to  parameters  passed  by  calling  routine,  and  starts  display.  Called  by  CONTL. 

13.  EFIX  - uses  INTEC  to  determine  coordinates  of  fixation  point  upon  canopy 
surface  for  pilot’s  nominal  eye  position  and  viewing  position.  Called  by  PERP. 

14.  PLAC  - converts  coordinates  of  vertex  in  object  space  to  their  equivalent  values  in 
perspective  space.  The  origin  of  this  space  is  the  pilot’s  nominal  eye  position  and  the  y-axis  is 
along  his  nominal  viewing  direction.  Called  by  PERP  and  TRSDC. 

15.  TRSDC  - converts  coordinates  of  surface  vertices  from  object  space  to  equivalent 
values  in  display  space.  Called  by  PERP. 

16.  INTEP  - removes  portions  of  figure  edges  which  lie  behind  viewer.  Called  by 

TRSDC. 

17.  TEST  - determines  whether  surface  vertex  is  outside  of  display  space  viewbox. 
Called  by  TRSDC. 

18.  CLIPL  - computes  apparent  intersection  points  of  test  edge  line  with  masking 
surfaces  and  generates  visible  portion  on  display.  Called  by  PERP. 

19.  CHKS  - computes  portion  of  test  edge  contained  within  displav  spave  viewbox. 
Called  by  CLIPL. 

20.  ORDLN  - orders  apparent  intersection  points  of  test  edge  with  edges  of  masking 
surface.  Called  by  CLIPL. 

21.  CKLIN  - determines  visible  portions  of  test  edge  remaining  alter  surface  masking. 
Called  by  CLIPL. 

21.  GIN,  GBEG,  GPUT,  GEON  and  GSTART  are  Fortran  graphics  routines  SEEK, 
GETDEV  and  GEVT  are  system  file  control  routines.  See  reference  6 for  details. 
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/ Joe(^'C^^ ) 

/CREiTElFINKSt-) 

/FORT 

C R4Y  TK4CING  FCR  PILCT 

C 3 ClRENSICN»t  FLAT  PLATE  CANOPY  WITH  OPSTRUCTICAS 
SUeROLTI^E  CtNTL 
CORPON/FCBG/LFG 
COBPCN/GABEA  / IGF IL ( ICCCO  » 

COBPON/C4W/NT,NP,NC,NB,NV(ICO),PXV(1CO,5»,PYV(ICC,5),P^Y(ICO,5) 

COFPON/NCRP/AXN(lCCt,AYN(lCC(,AZN(lCCI 

COPPON/L  INE/AStBS,GS,XS,YS,ZS 

COPPON/PILCT/XC,YC,ZC 

CIPENSION  AN(2) 

CATA  AN/2FYSt2HNO/ 

CATA  AIl,AI2/2.,2./ 

CATA  PIA3.1A159/ 

CALL  GETCEV(LLN, 'F/NKSH' ,1C) 

CALL  GETrEVUFG, 'FORGRA*  ,-l  I 
WRITE(1,S9P» 

998  FCRHAT( 2X, 'CANCPY  CATA  ON  FILFt  (YS  OR  NO)') 

REAC(  l.lCCl  )A 

IF  ( A.tC.  AN ( 1 ) ) GO  TC  1 
CALL  REACV 
CALL  NOR,«L 
GO  TO  2 

1 CONTINUE 
CALL  REACF 

2 CONTINUE 

WR  ITEI  I,  ICCC  ) 

ICCC  FCRPATI 2X , *Pk  INTOLT  CANOPY  CATAt  (YS  OR  NO)') 

BEACH, ICCl  )A 
ICCl  FORPATI 1A2 ) 

IF(  A .£(3.  AN  ( 1 ) ) CALL  TELTY 
WR  ITEI I,  :CC2  ) 

IOC2  F0RPAr(2X, 'CISPLAY  CANOPY  C CNF  I GUR  A T ! ON  t (YS  OR  NO') 

BEAU  l.lCCl  )A 
IF(A.EC.aN(2)  ) go  TC  3 
CALL  CRWCN 
WRITE!  1,  1CC3) 

10C3  FOBPAT(2X, 'PARC  CCPYS  (YS  CR  NO') 

REAC(  1,  ICCI  )A 
IF(A.EC.AN(2)  ) GO  TC  3 
WR  ITE ( 1,  ICCA  ) 

lOCA  FOBPAT ( 2X  , 'PR  INTER  , TURN  CN  ' ) 

PAUSE 

CALL  Gsnm 

3 CONTINUE 

WR  ITE 1 1,  icc: ) 

10C5  FCRPAT(2X,'CISPLAY  PERPFCTIVET  (YS  OP  NO)') 

REACI  l.lCCl  )A 

IF ( A .EG.AN ( 2 ) ) GO  TC  A 

CALL  PFRP(C) 

WR  ITE ( 1 , iCC J ) 

REACIltlCCnA 

IF!  A .EC.  AN ( 1 ) ) CALL  GS31  ( 1 ) 

A CONTINUE 

CALL  SEcK(LLN,C) 

RE  AC ( LUN  )NVW 
WR  I TE  ( 1,999  ) 

999  FCRPATI 2X, 'CCPPUT;  PRIMARY  RFFLFCTICN  PCINTSt  (YS  OR  NO)') 
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RE*r(l,lCCl)fl 
IF(fl.fcC.AN(2)  ) GO  TC  3CO 
JPC  = C 

V.«  HF  ( 1 ,SS6  ) 

■JSF  FOPRAI  ( 2X,  < Pk  INTOL  I C CMP  0 1 6 T 1 C K S F (VS  OR  NO't 
KEAr<l,lCCl)A 
IF ( A .EC.  1N( n ) 1PC  = 1 
NVLL=0 

CALL  SF':k  (LLA  ,\VW  I 
UftlTF (LLN  )NVLL 

C INCEX  PILOT  \/lEV.lNG  CIRECTION 

C fll  ELEVATION,  AN^.LE  FROM  I AXIS  TCFARl)  X AXIS  lA  X Z PLANE 
C AZ  AZIMLTH,  ANGLE  FROM  X Z OLANE  TChARC  V AXIS 

C Z AXIS  TCV«ARC  UPnARC,  X AXIS  TCWART  LEFT  FACING  FRONT  OF  CANOPY  AND  V AXIS 
C TCWARC  PACK  CF  CANCPY  ALONG  LCNGITLCINAL  AXIS  OF  AIRCRAFT 
A1=G. 

10  A2=AI2 

11  A2  = A2-U; 

IF( A2.GE.-9C.  ) GO  TC  15 
A1  = A1*AI  : 

IF ( A I .GT  . 18C  . I GO  TC  259 

RR ITE ( 1,SS5  ) A 1 
995  F0RNAT(2X,FlC.Pl 
GO  TO  10 
15  CCNTINUF 

air  = a;*p  I/IEl  . 

A2R  = A2*P  I/ieO. 

C CETERMINF  GIRECTICNAl  CCSINFS  CF  VIFWING  DIRECTION 
AS=CCS(A2R)*S1N(A1M) 

PS  = S IN( A2R  ) 

CS  = CCS<  A2M*CCS(  A IM  ) 

R = l. 

XS  = XC 
YS  = YC 

zs  = zo 

I PLN=C 

C CETERMINF  SURFACE  IMFRCEPT 
IT  CCNTINur 

FO  20  IS^l.NP 
ISK  = C 

CALL  INTECI  ISK, 1$,XR,YR,ZR) 

1F(  ISK.GT.C  ) GC  TC  25 
20  CONTINUF 
GC  TC  11 

25  CONTINUE 

IF(  IS.GT  .NP  ) GC  TC  26 
CO  TC  2C 

26  CCNTINUF 
IP(IS.LE.NC)  GC  TC  11 
CALL  CCRP(ANC,IS,Rr,TT) 

T=TT*R 

P=PT*R 

IF(R.LT..CCCC1)  GC  TC  11 
IPLN  = IPU".il 
A I*-AS 
F I=-PS 
C I =-CS 

PX*AS*AXr(lS)FeS*AYN(!S)+CSAAZMIS) 

AR--AS  + 2 .♦RXAAXM  I S I 
PR=-PS*  2 .ARX* AYN (IS) 
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CR*-CS  + 2.»Rx*ftZN(  IS) 

1F( IPUN.tC.2)  GC  TC  3C 

IPP  MS 

*SP=AS 

PSP  = es 

CSPMS 

XPslCP 
YP  = YR 
/ P = / « 

^ IPM  I 
P IP  = P I 
CIP=CI 
«RP=flR 
FRP=PR 
CPP=CR 
dS=-«K 
PS=-PH 
CS=-CP 
XS  = XR 
YS  = YR 
IS  = r\i 
CC  TC  17 
3C  CONT  IFU' 

N VLL  =MVLl < 1 

WRITf(LLi\)«l,A2,ASP,PSP,CSP,IFP,XP,YP,2P,AIP,PIP,CIP,ARP,PRP,CRP, 
CRfIS,XR,YR,ZR,fll,ei,CI,4R,BR,CR,T 
IF ( IPC  .FG  .C  ) GC  TC  11 
V.RITEn,PCC)fll,'\2,aSP,°SP,CSP 

WRITFI  5f  CniPP.XP.YP./P.MP.PlP.CIP.RRP.HRP.CRP.P 
V»RITF(3,'C1)IS,XR,YR,2R,AI,PI,CI,AR,PR,''.R,7 
ecc  FCPPRT(2<,2(F7.2,2X),3(F7.^,?X)) 

PCI  FCR»'ATI2x,I3,3ir6.2»2X),7(F7.R,2X)) 

CC  TC  11 
2SS  CCNTIMJF 

call  SSi'k  I LLN  .NVV.  ) 

k.R  ITtILL',  )NVLL 
ICC  CCNTINUT 

WR  ITF I .CC6  ) 

ICC6  FORMAT  ( 2X,  • PtRPcCT  I V£  ' /2X  , 'RAY  F'jTRFNCE  PC  I N T S ' / 2 X , • R I GP  T HAND  SIC 
CF  CCCKPI  ! ( YS  CR  NC  ) * ) 

R F A r ( 1 , 1 C C 1 ) A 

IF  ( '.  .tC  . 1M  2 ) ) GO  TC  AC 

CALI  PFRP(l) 

)« R I T F ( 1 , : C C 3 ) 

PFAni.ICCl  )7 

IF  ( • .EC  . 'M  1 ) ) CALL  GS  31 ( 1 ) 

«C  CCATINLE 

RRI’EH,  .CC7) 

10C7  FCRVAT(2X, ‘PERPRCTIYF '/2X,  'PRIMARY  R F F LFC T I 0 N S t • ) 

REAr(l,lCCl)« 

IF  ( A .ff.  . 'M  2 ) ) CO  TC  A1 
CALI  PF-iPt2) 

R R I T E ( 1 , 1 C C 3 ) 

REACI  l.UCl  )A 

IF(  A .iC  . AM  I ) ) CALL  GS  31  ( 1 ) 

«1  CCNTINLF 

CALL  CFV r (LLN ,C,C  ) 

PETLRN 

FNC 

/FCPT 
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00000'S*  ooo 


IN  SLSF'KF  vERTlCFS 
SLEKCLTl'.E  RE&Cv/ 

ccnvcn/C'n/m,np,nc,np,\'v(1CC>,pxv(:co,‘--),pyviicc,5),p/v(ico,si 
CCPVCN/V  BT/XV(:0(),YV(lCC),ZV(:CUi),«VR(’.CC,E) 
CCf'PCN/\(RP//'XMlCC),4YN(lCC),&Z‘l(lCCI 
CCPPCN/P llCT/XC,YC,ZC 
P E 4 r ( 2 , 1 C C C ( 

ICCO  ECP«'4T(  I 

PE4L'(2,  ICCl  )M  ,’'IP  ,NC  ,NP 
lOCl  FCP“4T ( /P ( 2X , 12  ) ) 

PE#r ( 2,  ICCC  ) 

PE4ri  2,  1CC2  MNV  ( I I , I = l ,'JPI 
1CC2  FCPYflTI 1F(2X, 12) ) 

PE4E  ( 2,  UCC  ) 

CC  1C  .1=1, E 

PE4r(2,lCC2)(NVP( I,J),I=L,KP) 
iO  CONTINU'- 

PE4r<  2,  UCC  ) 

PE4C  (2,  ICn  ) ( XV  ( I ) , 1 =1  .NT) 

10C3  FOPP4T(P(2X,F7.R)  ) 

PE4r(2,UCC) 

RE4C( 2, ICC3  M YV( U , 1 = 1 ,NT  ) 

RE4C  ( 2,  UCC  ) 

BE4r(2,lCC?)(ZV(  I ),1  = 1,NT) 

BE#r ( 2, ICCC ) 

RE4C(2,1CC3)XC,Y0,ZC 

RETLRN 

ENC 


EST4eLl.SF  SUPF4Cr  NCPN4L  FOR  EACH  PLATE  SURFACE 
SURFACE  NCRP'LS  CIRECTEC  TOWARD  CCCKPIT  INTFRIOR 
SUERCLTINF  NCRPL 

CCPwON/CAN/NT.NP.NC.ND.NVdCCl.PXVllCO.Sl.PYVdCC.Sl.PZVIlOO.S) 

CCf'PCN/VCRT/XV<lCC),YVdCC).ZVdC0).NVRdCC,5) 

CCPNC\/NCRR/AXM:CO),AYNdCC),AZNdCO) 

COPPCN/P  ILCT/XC.YC.ZC 
NVW  = C 

CALL  GFTI'EV  ( LLN  , • F INKSH' , 1C  ) 

CALL  SEFK(LLA.C) 

WR  ITE (LLN  )NVW 
WRITE (LLN  )NP ,NC ,NP 
CC  :C  1 = 1, NP 
NK  = KV(  1 ) 

WR  ITf (LLN  >NK 
CC  = K=1,NK 
KV  = NVR(  I ,K  ) 

P>V(1,K)=XV(KV) 

PYV ( I ,K ) = YV (KV  I 
PZV(  I,K)=ZVtKV) 

WRITEILLN)PXV(l,K),PYVd,K),PZVd,K) 

5 CCNTINLF 
K = 2 

7 A1=PXV(  I,X  )-PXV<  1,1) 

A2  = PXV(  I ,K+  1 )-PXV ( I , 1 ) 

P1=PYV( I ,K )-PYV ( I , 1 ) 
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O O O >V  o o o 


F2=PYV(I,K4l)-PYV(l,l) 

C1  = P2VM  ,K  )-P2V  ( I , I I 
C2-‘'2V(I,K«1)-P2V(1,1) 
P!=SURKJl*»2^Rl*Y2>Cl»*2t 
P2=SCBT( Y2*»2«P2**2*C2**2) 
/)-(/>I*/S2*ei»P2<Ci*C2)/(Pl*P2) 
IF  ( ■'PS(  'U  .1  I . 1.  ) GC  TC  9 

K i«  < 1 

IF(K.tG.\K)  (C  TO  iC 
CC  TC  / 

9 iN  = !»CCS(AI 

R = 1 ./ ( p: »P2»S  IN ( AN  ) ) 
iXN ( I »=- ( P 1*C2-C  l*H2  )*R 
AYN ( I) ( A 1*C2-C 1*A2 ) *9 
AZM  1 1 =-  I Al^P2-A2*Hl  ) *R 
WRlTf(U.N)AXMn,AYN(I),AZN(n 
1C  CCNTINOF 

1.R  ITE  (LLN  )XC  , YC  ,Zt 
NVW=N£XR  C(LLN) 

CALL  SEFK(LLN,C) 

hR  I TF  ( LLN  INVVi 

CALI  C£V I (LLN tC,0 I 

RETLRN 

ENC 

/FORT 


WEAC  FUE  FCR  CANOPY  CATA 
SLERCLTINE  REACF 

CCRRCN/CAN/NT,N?,NCtNP,NV(lC0),PXV(in0,''i),PYV(lCC,5),PZV(100,5l 
CCRVC'N/N'CRP/AXN  ( ICC  ) , AYN  ( ICC  ) ,AZN  ( ICO) 

CCRRCN/PlLCT/XC.YCtZC 

CALL  GFrCEVULN, 'FINKSH' ,1C) 

CALL  SE>^K(LLN,C) 

REATILUN  )NVk« 

REACtLljN  )NP,NC,NP 
CC  :C  1^1, NP 
PEACtUJN  INK 
NV ( I ) =NK 
CC  E K=1,NK 

REAC(LL)N)PXY(I,K),PYV(I,K),PZV(1,K) 

5 CONTINUE 

REAC(LLN)AXN(I),AYN(I),AZMI) 

1C  CCNT  INU' 

PE  AC (LLN  )XC  , YC,/C 
call  CEVT (LLN ,C  ,C  ) 

RETLRN 

ENC 

FORT 


PRINTCLT  CF  CANCPY  DATA  FCR  REVIFK 
SUP9CLTINE  TELTY 

CCRR0N/C'N/NT,NP,NC,NP,NV(1CC),PXV(100,S),PYV(1CC,5),PZV( 100,5  ) 

CCRRCN/NCRR/AXN(ICC),AYN(1CC),AZN(ICCI 

CCRRON/PILCT/XC.YL.ZC 

V.RITt(  USSR) 

998  FCRRAT( 2X  , *PK  INTER  , TL«N  CN  ' ) 

PALSE 

V.R  ITr  ( 3,  ICCC  )NT  ,NP,NC  ,NB 


ooon>.nr>  o orio 


ICCC  FCPvaI1A(2X,12)) 

CO  1C  I=:,NP 
KN  = NV  ( I) 

CC  5 K=1,KA 

WRITt(3,:CCnPXV(l,K),PYV(I,K),P2V(l,K) 

1001  FCRPAT(2X,2(F6.2,2X)) 

5 CCMINCF 

WRlTEn,:CC2)AXMn,AYN(l),A2K(l) 

1002  PCRPAT( 3( 2X,F  IC.^  ) ) 

1C  CCMINOF 

WR  ITE  < 3t  UC2  )XC,YC  ,2C 

REICRN 

ENC 

/FtPT 


CETERPINES  IF  RAY  STRIKES  CONVEX  SLREACE 
SL'ERCcriNE  INTECI  ISK,IS,XR,YR,2R) 

CCPP0N/C5N/M,NP,NC,NP,  .NVdCCl.PXVdOO.ll.PYYdCC.SI.P/VIlOO.S) 
CCI-PCN/NCRR/AXM.CO.AYNdCO.AZFJdtC') 

COPF'CN/L  INE/AS,PS,CS,XS,VS,2S 

CK  = AXMIS)»ASFAYMlSI*fTS  + AZMlS)*CS 

IFICK.GC.C.)  RETURN 

HAY  STRIKES  SURFACE  IN  CLTWAKC  CIRECTICN 
1=1 

IF(XS.FQ.PXYdStn.AND.YS.EC.FYVdS,n.ANC.ZS.EC.FZVlIS,I)d  = I + l 
S=(AXMIS)*(PXV(IS,I)-XS)+AYM15)«(PYVdS,II-YS)+AZMISI*(PZV(IS,I 
C l-ZS  ) IZCK 
XR=AS*S*XS 
YR=aSyS«YS 
ZR=CS«S«/S 
IN  = NV(  IS  ) 

CC  :c  i=:,iN 
ic=id 

IFII.EC.IM  fC  = i 
A1 =PXV ( I S , I )-XS 
e i = PYv( d , I )-YS 
C1  = PZV(  I S,  1 )-ZS 
A2  = PXVdS,IC)-X$ 

P2=PYV(  IS,  IC  )-YS 
C2=PZV(  IS,  IC  )-ZS 
p 1 = XR-XS 
P2=YR-YS 
P3=ZR-ZS 

C = Pl*(",l«'C2-H2*Cl)-P2*(Al«C2-Cl*A2dP3*(Al*P2-Pl»A2) 

IFIt.LT.C.)  RETURN 

KAY  STRIKtS  SURFACE  ON  FNCUOSED  SIPE  OF  SURFACF  ECGF 
1C  CCNT  INUE 

RAY  STRIKES  ENCUCSEC  SURFACE 
ISK  = 1 
RETURN 
ENC 

CRT 


CCKPUTES  INCICENCE  AfiGUE,  REFLECTArCE,  and  TRANSMTTANCt 
natural  UlGFT,  ACCITICN  CF  PCLAPIZAIION  CCYPCNCMS  IGNORfC 
SUPPCLTINF  CCPP(ANG,IS,RT,TT) 
CCPPCN/N^RP/AXMlCCdAYNdCCdAZN'dCCI 
CCPRCN/UNF/AS,PS,CS,XS,YS,ZS 
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C XN,  INCEX  Of  REFKACUCN,  TX,  IMERN6L  TRftNSNITTflNCE 
C4T/S  XN,  rx/1.5,  .9i/ 

*=-aS*aXN(lS)-BS*AYN(IS)-CS*AZN(IS) 

IF (ABSt « t .CT  .1.  I A = l. 

ANG  = ACOS( A ) 

ANrp=  as  IK ( S IN ( ANG ) /XN ) 

CA=CCS(ANG) 

SA  = S INI ANG  ) 

S1=SCRTIXN**2-SA**2) 

R0=-( ( (CA-Sl  )/(CA+Sl)  )**2*{ (CA*(XN**2)-S1)/(CA*IXK**2)+S1>  )**2)/2. 
T0  = ( l.-RO  )*CA/CCS( ANGP  ) 

CA^CCSI ANGP  ) 

SA=SIN(4KGP) 

Sl  = SCRT(  XN**2-SAX'*2) 

RI  = (((Ca-Sll/(CA-fSl)  (CA#(XN**2»-S1)/(CA*(XK**2)+S1)  )**2)/2. 

TI  = ( l.-R  I )*Ca/CCSI4NG) 

TT  = TO’*TMTX/(l.-(RI*TX)'**2) 

RT=RO+R I*TT 

RETURN 

ENG 

/FORT 

C 

C 

C CCRPITES  ARCSIN 

FUNCTION  ASIMXI 

Y = X 

AX  = ABS( Y ) 

IF(AX.GE.l.C)  GC  TC  A 
AC^ATAN'IY/SCRTII.-Y^Y)) 

AS  IK  = AC 
RETURN 

A IF ( AX.GE  . l.CCOl  ) GO  TO  10 
AS  IN  = 1.5  IC7S 
IF(Y.LT.C.)ASIN=-ASIN 
RETURN 

10  hRITEIl.lDX 
STCP 

11  FCRPATISFERRCR  ♦ ,23HARCSIN  ARGUMENT  .GT.l.  ,6HAKC  = ,f:6.8I 
ENC 

/FORT 

C 

C 

C PRL  ARCCOS  ROUTINE 
C REF.  GLENN  BECK 

FUNCTION  ACCSIX) 

Y = X 

AX  = aeS( Y I 

IFIAX.GE.l.CIGCTO  A 

IFIY.cO.C.  IGCTC  3 

AC  = ATAN( SCRT ( l.C-Y*Y ) /Y) 

1 IFIY.LT.C. IGCTC  2 
ACCS=AC 

RETURN 

2 ACCS=AC» 3.lRi5926 
RETURN 

3 ACCS=1.5<C7St3 
RETURN 

A IFIAX.GE.l.CCCClIGCTC  10 
AC»C  . 

COTC  1 


non 


10  WR  ITE(  1.  1 llx 
STOP 

11  FORRATI  PEERRCR  ♦ ,23HflRCCCS  fiRGLRtNT  .GT.l.  ,6H/^^C  = ,ei6.8l 
ENC 

/FORT 


CRflhS  GRAPHIC  PICTCRF  CF  CANOPY  IN  3 FOLD  LAYOUT 
SUERCLTI‘E  CRhCN 
COPPCN/FCRG/LFG 
CCPPCN/CAREA/ICFILdCCCO) 

COPPON/CAN/M.NP.NC.NP.NVdCO.PXVdCO.SI.PYvaCC.SI.PZVdOO.Sl 
CCPPCN/NCRP/AXN(  ICC  ) lAYNdCC  ) ,A/NdCO) 

CCPPCN/PIICT/XC,YC,ZC 
CATA  SX/C.C235/ 

WRITE!  l.SSe  ) 

q<)8  FCRNAT  ( 2X  , 'CCNSCLF  NO  1,  TURN  CN»| 

PAUSE 

CALL  GIMFCCC) 

IX=(YC-5C.)+SX 
IY=(ZC-eC.  )*EX 
CALL  GPEG!  1,  IX,  lY  I 
CALL  EPA-iK 
IY=( XC  + RC  . )*SX 
CALL  GCPY(2,;,IX,  lY) 

1X=(ZC  + 8C.  )*SX 
CALL  GPEGO,  IX,  lY) 

CALL  ETIC 
NE  = 3 
NS=N6+1 
CC  ;c  I=NS,NP 
C FNTITY  IS  CANCPY  SURFACE — FENCE  CR  TRANSPARANT 
CS  = -AXN(  I ) 

CF=-AYN(  I ( 

CT  = AZN(  I ) 

IX  = ( PYV(  1 , 1 )-5C.  ) «SX 

IF  I ( I .LE  .NC  .ANC  .OS.LT.O.  ) .CR  . d .GT.NC.  AA'P.CS.GT  ,C.  ) ) CO  TC  2 
C SURFACE  FACES  VIEWER  FRCP  SIDE  VIEW 
IY=(PZV(  1,1  )-6C.  >*SX 
NE=NE+1 

CALL  GPcClNE,  IX,  lY  I 

IFCl  .GT.NP.ANC.I.LE.NOCALL  GPLT(3,133,1,2> 

CALL  L INS (I  ) 

2 CONTINUE 

iY  = ( pxv!  I , n + sc . I *sx 

IF ( ( I .LF .NC .ANC.CT.LT.O. ) .OR. (I .OT.KC. AND.CT .GT.C.  ) ) GO  TC  3 
C SURFACE  FACES  VIEWER  FRCP  TOP  VIEW 
NE=NE*I 

CALI  GHEG (NE  , I X , I Y ) 

IFd.GT.NP.AI.C.I.LE.NOCALL  GPLT(3,130,1,2) 

CALL  LINT!  I) 

3 CCNTINUr 

IF ! ! I .LE.NC .ANC .of .LT .0.  ) .CR.  ! I .GT.NC. ANn.CF .GT.C.  ) ) GO  TC  10 
C SURFACE  FACES  VIEWER  FRCP  FRONT  VIEW 
IX  = !PZV!  1,1  )a5C.  )*SX 
NE=NEA1 

CALL  GPFC|NE,IX,1YI 

IF! I .GT.NP.ANC. I .LF.NC  TCALL  CFLT  I 3 , 1 30 , 1 ,2 ) 

CALL  L INF  ! I ) 

10  CONTINUE 
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■ 


CC  20  1 = 1, Ne 
CALL  GcCMl) 

20  CCNTINLiF 
CALI  CSTART 
RETLPfN 
ENC 

/FORT 

C 

C 

C mark  eve  PCSniCN  IN  CANOPY 
SLERCLTINE  EMARK 
CALL  GPUM6,EC,-lC,-5) 

CALI  GPLT n t 7C,C,  1C  I 
CALL  GPLI(8,EO,10,-5I 
CALL  GPLlT{<5,7C,-7,-lC) 

CALL  GPLT(1C,5C,0,2CI 

RETLRN 

ENC 

/FORT 

C 

C 

SUERCLTINE  ETIC 
CALL  GPL7 (t,7C,C,-5) 

CALL  GPL  7 (7, EC, 0,10) 

CALL  GPL M e, 70, -5,-5  ) 

CALL  GPL7(5, 50,10,0 

RETLRN 

ENC 

/FORT 

C 

C 

SUERCLTINE  L1NS(  I ) 

CCMMCN/CAN/NT,NP,NC,NB,NV(1CCI,PXV(100,5),PYY(1CC,5),PZV(100,5) 
CATA  SX/6.C225/ 

CALL  GPUTO,  13C,2,C) 

NK=NV( I ) 

CC  10  K=1,NK 
K 1 = K + 1 

IF(K.EC.NK)  Kl=l 
1X=(PYV(  I,K1)-PYV(I,K))*SX 
IY=(PZV(1,K1)-PZV(1,K))*SX 
CALL  GPLT(K*5,53,1X,IY) 

10  CONTINLF 
RETLRN 
ENC 

/FORT 

C 

C 

SUERCLTINE  LINF(I) 

CCRMCNZCAN/NT,NP,NC,Ne,NV(lCCI,PXV(lCO,5),PYV(KC,5),PZV(100,5) 
CATA  SX/6.C235/ 

CALL  GPLT(3,13C,2,C) 

NK=NV(  I ) 

CC  10  K=1,NK 
K1=K41 

IF(K.EO.NK)  Kl=l 

ix  = (pzv(  i,K 1 )-pzv( i,K n*sx 

IY=(PXV(1,K1)-PXV{I,K))*SX 
CALL  GPLT(K45,53,IX,IYI 
10  CONTINUE 


J 


o o o n >%  no 


BETL«N 

EKC 


/FOPT 


SUFPCLT  r.e  L INT ( n 

CCPPCK/CJN/M,NP,NC,NP,NV(lC0),PXV(100,‘j),PYV(lCC,5»,PZV(100,5» 
CAT/  SX/f.Cii5/ 

CALC  GPLT(3,13C,^,CI 
NK=NV ( I ) 

C C : C K = . , N K 
K1=K*1 

IF(K.EC.^K)  Kl=l 
1X=  ( PYV(  I ,K  1 )-PYV(  1 ,K  n*sx 
IY=(PXV(  I,Kn-PXV(I,KI)»SX 
CALL  GPL T (K«E ,53 , I X , I Y » 

1C  CCMINLE 
PETLPN 
E\C 


LBAWS  PFpPrCIlVF  CF  CANCPY  ANC  EKTBY  AKC  REFLECTIOS  FOI^TS  PILOT  EYE 
PCS  IT  IGF 

SLEPCLTINE  PFPP(ILS) 

CCPKCN/FCRG/LFG 
CCKPCN/GAREA/ICFILIUCCO) 

CCPPCN/P'R/NC,NS(lCCI,PXS(UC,lC),PYS(no,10»,PZS(lGO,lOI 
CATA  AX, PX,CX/1C . , 5.  , 512.  / 

CALL  TRSCC 

CALL  CFTCEVILLN, 'FISKSH' ,1C) 

CALL  GIN(ICCCC) 

CALL  EE  IX (XX , V Y ,ZZ  ) 

CALL  PLACIYY,ZZ,PY,PZ  ) 

IXE=XX*AX*CX/(PY*exl+CX 
IYE=PZ*AX«CX/(PY*eX)«CX 
CALL  GPFC(l,IXE,IYt» 

CALL  ETIC 
CALL  GFC*-(1» 

NE  = : 

IKMC 

cc  1C  i=:,NC 
NK  = 5S(  I) 

IF(lK.r.T.t)  CALL  GBFG  ( NF  , I XF  , I YE  ) 

IK  = F 

CC  ' K=:,NK 

K 1 = K 4 1 

1F(K.FC.‘>K)  Kl=l 
C CEECK  ECGES  FCR  F ICCEA  LINES 
CALL  CLIPLI  IK.I.K.Kl) 

5 CCNT  INUF 

IF ( IK  .CC  .e  ) CC  TO  10 
CALL  GEC^  (NE  I 

NE=NE4l 

1C  CCNTINLF 

IFI  lES.EC.CI  CC  rc  21 
CALL  GHEC(NE41,IXE,1YF) 

CALL  GPLI(5,17eC,C,.FALSF.I 

IFI IcS.EC.2 ) CALL  GPL T I 5 , 1 75C ,C ,0 ) 

CALL  SFFKILLN.CI 


■it  rr  * 


o o u 


REiC (LLN  INVRC 
CALL  SEEK (LLN »NVwC) 

RE«C(LLIN  )^VR 
1K  = 6 

cc  ic  i=:,NVR 

RE<C(LUN»/llV,A2V,flSP.RSP,CSP,IP,XP,yP,Z0,»IP,eiF,ClF,«RP,eRP,CKP 
CR.ISiXR,VR,ZR,iI,Hl,CI,4R,BR,CR,T 
IF(  lES.FC.l  I GC  TC  3C 
GC  TC  AC 
2C  CCNTIRUF 

CALL  GCCN(NE<1) 

21  CALL  GSTART 

CALL  CEVT ( LLN  ,C ,0  I 
PETLRN 

C FNTRENCF  PCIMS  FCR  EXTERNAL  RATS 
3C  CCNT  INIJF 

CALL  PLAC(VR,ZRfPY,PZ  I 
IFJPY.LE.C.)  GC  TC  2C 
XX  = XR*AX/ (PY«BX  ) 

22=PZ*AX/(PY»BX» 

IF ( ABS( XX  1 ,G I . 1 . .CW  .ABS( ZZ  ) .GT.  1 . ) GC  TO  20 

1X=XX*CX4CX 

IY=ZZ*CX«CX 

CALL  GPUT( IK, A3, IX,  lY) 

IK  = IK4  I 
GC  TC  2C 

C PRIRAHY  REFLFCTICN  POINTS 
AC  CCNTINLIF 

CALL  PLAC(YP,ZP,PY,PZ) 

IF ( PY  .L?  .C  . ) GC  TC  2C 
XX  = XP*AX/(PY*BX  ) 

ZZ*PZAAX/ (PYABX  ) 

IF<ABSIXX).GT.l..CR.ABS(ZZ).GT.l.)  CC  TO  20 

IX=XX*CX«CX 

IY=ZZACX*CX 

CALL  GPL T ( IK, ICC, IX-A,C) 

CALL  GPL  I(1K+ 1,  lie, IY-A,C) 

1C=-ALCG  1C ( T I 
IC=  IC4176 

CALL  GPLTI IK42,9C, 1C,C) 

IK  = IK43 
CC  TC  2C 
ENC 

/FORT 


CCRPLTES  EYE  FIXATION  CENTER  FCR  FFRPFCTIVE  DRAWING 
SLBRCLTINE  £F  IX  I XP,YR,ZR  ) 

CORRCN/CAN/NT  ,NP,NC,NB  ,NV( ICC  ) ,PXV(  1C0,S)  ,PYV ( 1CC,5  I ,PZV(  100,5  I 

CCRRCN/P1LCT/XC,YC,ZC 

CCRRCN/LINE/AS,eS,CS,XS,YS,ZS 

CATA  AN,PI/5.,3.1A15'5/ 

ANN=AN*PI/iec. 

XS  = xc 
YS  = YC 
ZS=/C 
AS'C  . 

PS  = -CCS(  'NN  I 
CS  = S INI  a,'  N ) 

cc  2C  IS--1,NP 


o n n r>  N*  o r»  o n o 


1SK  = C 

C#LL  INTeC( ISK,IS,XR,VR,ZR) 
IFdSK.GT.Cl  RETURN 
20  CONTINUE 
RETURN 

EMC 

/rOMT 


CONVERTS  OBJECT  SPACE  CCCRDINATES  IMG  EVE  SPACE  CCCRCINATES 

PILOT  EVE  CIRECTED  IN  Y-Z  PLANE  CF  AIRCRAFT  CCORCINATES  AND  ‘>-CEGREFS 

APOVE  NEGATIVE  V-AXIS 

SUEROUTINE  PLAC(PYV,PZV,PV,PZ  I 

COPPCN/PILCT/XC,YC,ZC 

DATA  AN/S./ 

CATA  PI/3.1A15S/ 

AN1  = AN*P  I/ieC. 

PV  = -(PVV-YC  )*CCSIANn-(PZV-ZG)*SIN(AM) 

PZ=-IPYV-YC  IASINI  ANl  H(PZV-ZC  )♦COStA^U 

RETURN 

ENC 


CONVERTS  OBJECT  COORDINATES  TC  DISPLAY  CCORDINATtS  , REFCVES  LINES  BEHIND 
VIEWER.  SURFACES  OUTSIDE  VIEWING  BCX,  SURFACES  NOT  F/CING  VIEWER 
SUBROUTINE  TRSCC 

C0f<H0N/CAN/NT,NP,NC,N6,NV(lC0)  .PXVdOO.S)  ,PYV(1CC,5)  ,PZVl  100,5) 
CCRWON/NCRR/AXN(  ICO.AYNdCO.AZNdCC) 

CCPPCN/PILCTAXC.YC.ZC 
COPHON/L  INE/X1,Y1,Z1,X2,Y2,ZZ 

COPPON/PFR/ND.NSI  lCC),PXSdCC,lC),PYSdOO,lO),PZSdOO,10) 

CIPFNSICN  XS(lCC,5),YSdCC,E),ZSdC0,5) 

CATA  A,B/1C.,5./ 

NST  = NB*  I 
NC=  1 

CO  10  I=NST,NP 

C»AXMI)A(PXV(I,l)-XC)+AYMI)A(PYVd,l)-Y0l*4ZNd)A(PZVd,  1 )-Z0) 
IFIC.GT.C.)  GC  TC  IC 
C SURFACE  FACES  VIEWER 
KN*NV(  I ) 

CC  2 K^l.KN 
XSINC.K  ) =PXV ( I ,K ) 

CALL  PLAC(PYV(I,K),PZVd,K),YS(NP,K),ZS(NC,Kd 
C CCNVERTEC  TC  VIEWER  CCCRDINATES 
2 CONTINUF 
KC»C 
IXP-C 
IXN  = C 
IZP  = C 
I Z = 0 

CO  5 K=1,KN 
K1*K41 

IF  (K  .EC  .KN  ) K 1=  1 

IF ( YSINC  ,K  ) .LT. C . . AND.  YS (NC  ,K1  ) .LT.C.  ) GO  TO  5 
IF ( YSINC ,K  I .L T.C.  .CR.  YSINC.Kl  ) . L T.C.  ) GO  TO  3 
X1  = XS(NC,K  ) 

Y1*YS(NC,K  ) 

Z1=ZS(ND,K) 

C=il 
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r 


GC  TO 

3 COMIMJE 

C«LL  INTlP(XS(N^,K),YStND|KI,^$(Nn,KI,XS(^C,^<n,YS(^C,Kll,/S^^•C,•tl 
C )»C  ) 

C SURFACE  ECGE  EXTENCING  REMND  VIEWER  TRIRMEn 
t,  KC  = Kr«  i 

PXS(NC,KC)=X1*A/(V1*R) 

PYS(NC,Xr  )=-l./Yl 
PZS(NT,KC  ) = Z1*A/(Y1*P) 

CALL  TFST(PXS(NC,Kt)),PZS(^D,KCI,IXP,IXN,l/P,U^) 

C PCSmvf  PO(M  Cfi  CROSSING  LINE  FRO  BEHNO 
IFK.GT.C.)  CC  TO  S 
KC»Xf*l 

PXS (NC.Kt  l = x:*A/ ( V2*B ) 

PYSINC.Kt 

PZS(NC,KC)»Z2*A/(Y2*RI 

C LFAC  PCINT  FCP  LINE  CROSSING  FRCP  IN  FRCM 

CALL  TESHPXSINC,KC»,PZS(NC,Kn),IXP,IXN,IZP,IZN) 

C VFRTEX  CCNVFRTEC  TC  DISPLAY  VIEWDCX  COCBDINATFS 
5 CCNTINUF 

IFIKC.tC.C)  CC  TO  IC 
C SURFACE  NCT  REUNC  VIEWER 

iFi  ixp.ge.kc.cr.ixn.ge.kc.cr.izp.gf.kd.or.izn.ge.kc)  go  to  iO 

C SURFACE  PARTIALLY  C»  CCPPLETELY  wITUN  VIEW  nOX 
NSINC  »=K! 

NC=f  C* I 
IO  COM  INLE 
NC  = NC-  1 
RFTLRN 
ENC 

/FORT 

C 

c 

C TEST  SURFACE  VERTEX  f CR  POSITION  CLTSinf  VIFWPCX 
SURPOLT INE  TEST! X,Z  . I XP,  I XN  ,I ZP,I ZN ) 

: . ) 1 XP=  I XP*  I 

-l.)IXN=IXN*l 
1.)  1ZP=I/P«1 
-I  . ) IZN=  !/’,♦  1 


IFIX.CT, 

IFIX.LT, 

IFIZ.CT, 

I F ( Z .LT , 

RETURN 

ENC 


/FCRT 
C 
C 
C 


calculates  PTFRSFCTICN  PCINT  FCP  LINE  WITH  PLANE  N.OEPAL  TO  Y-AxIS 
SUErCLTI''E  IMEPIX1,y1,Z1,X2,YJ,Z?,R) 

CCPPCN/L  INE/»S,YS,/S,XF,Yf,ZF 
CAT  A YO/*./ 

R = SlPT((x:-X2)a»2MY1-Y2)»*2*IZ1-Z2»*»2I 
A = ( X2-X  1 > /R 

f = (y2-y:  »/p 
C = (/2-/:  l/R 
IF ( r .EG  .C  . ) return 
ECGE  NCT  PTRALLEI  TC  PLANE 
Y I » YC 

R I = ( Y I-Y  1 )/P 
X I = X 1 ♦ A*  I 

Zl-Zl«f*  I 
IFIT'.LT.C.)  LC  TC  3 
LINE  CPICP.ATES  PEEINC  VIEWER 
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n o r>  o 


XS  = X I 
YS  = V I 
/S  = /I 
XF  = *2 
YF=Y2 
2F  = /2 
RETLRN 
? CCMIMJK 

C lINF  CRIGUiaTES  IN  FRCNT  CF  VIEWER 
XS  = X1 
YS  = Y 1 
2S  = / 1 
XF  = X I 
Yf  = Y I 
ZF  = ; I 
RETLRN 
ENC 

/FORT 


TETERRINES  VISIPLE  PCRIICN  CF  ECGF  PY  RF^OVING  HICCFN  LINF 
SLL  CCNSTR-M'  T SLRFSCES  CONVEX 

SCPRCLTI'E  CLIPLI  IK.IV.KV.KVl) 

CCPPCN/lINE/AVfPV,CV,XV,YV,ZV 

CCRNCN/P'-fi/Nr  ,NS  ( ICC),PXS(lC0,lC»,PYS(ir)O,10),P2S(lCn,lC» 
COPPCN/CR/IW/NSC,RV,RCV(1C),RFV(10) 

CM/t  CX/‘  12  ./ 

XV  = PXS(  IV, KV  ) 

YV  = PYSt  IV, KV  ) 

ZV  = PZS{  I V,KV  ) 

RV  = SCRT((PXS(IV,KVl)-XV)**2-*(FYS(IV,KVn-YV)X*2+(FZb(IV,KVl)-ZVIw* 
C2  ) 

«V  = ( PXS(  IV ,K V 1 l-XV  l/R V 
ev  = IPYSI  I V ,K V 1 )-YV  I /RV 
CV=(PZS(  IV,KV1)-ZV)/RV 
RO  = r . 

RF=RV 

C/iLL  CHKS(-1.,C.,C.,R0,RF) 

CALL  CHKS(1.,C.,0.,RC,RF  I 
CiLL  CEKS(C.,C.,1.,RC,RF) 

CALL  CFKS(C.,C.,-1.,RC,RF) 

IF(RF.LT.PC)  RETLRN 
C eCGE  CFFCKEP  AGAINST  ViEWROX 
NSC  = 1 
RCVI  1 )=»C 
RFV(  1 )=RF 
CC  1C  I==1,NC 
IF  I IV. EC . 1)  rc  TO  1C 

C FINE  APPARANT  INTERSECTICN  POINT  CF  EOGFS 
ISC:=C 
KN  = NSII  ) 

CO  E K=1,KN 
Kl=K4l 

IF(K.FCi.KN)  Kl  = l 
XE  = PXS(  I ,K ) 

YE  = PYSt  I ,K  ) 

ZE  = PZS(  I ,K  ) 

RE  = SCRT((PXS(I,Kn-XE)**2*(PYS(I,Kn-rtl**2-*(PZSn,KI>-Zc)**2) 

AE= I PXS(  I ,K  1 l-XE  » /RE 
EE  = (PYS(  I,K1  )-YF)/RE 


. 1 


CE=(PiS(  I,K1  )-ZE)/RE 
C=<E*CV-AV*CE 
IF(C.EQ.C.)  GO  TO  5 
C ECGES  NOT  P4R«LLEL 

RIV=«CE*(XV-XE|-aE*(ZV-ZE)l/Q 
IF(*E.FC.C.  I GC  TC  2 
R IE=  ( XV-XE4/i\(*R  I V »/4E 
GC  TC  3 

2 CCNTIKUF 

PIE=(ZV-Z£'*CV*RIWI/CE 

3 CONTINUE 

YIE  = YE  + nF4R  IE 
YIV=YV*PV*RIY 

1F(Y1E.GI.(YIV-.CC1)  ) GO  TC  5 
C CONSTRAINT  ECGE  IN  FRCNT  OF  TEST  EOCF 

IF (R  lE.G  I .0  . .ANC  .R  I£  .LT.RE  ) GC  TO  7 
S CCNT  INUE 

IF(ISC.LT.2)  GC  TC  1C 
IF(RO.LT.C..ANC.RF.GT.RV)  RETURN 
C TEST  ECGE  NCT  EICDEN  PEEINO  CCNSTRAIMT  SURFACE 
IFIRC.GT.RV.CR.RF  .LT.C.  ) GC  TC  10 
C TEST  ECCE  PARTIALLY  BLOCKED  BY  CONSTRAINT  SURFACE 
call  CKLIN(RC,RF) 

IFINSC.EC.CI  RETURN 
GC  TC  1C 
7 CONTINUE 

C CONSTRAINT  ECGE  LOCATED 
ISC=  ISC+  * 

IF!  ISC. EG. 1 I R1=RIV 

IF(ISC.EG.2)  call  CRCLNIRl ,P I V.KO.RF  ) 

GC  TO  S 
1C  CONTINUE 

C VISIBLE  PCPTICN  CF  FOGE  REGAINS 
CC  2C  1=1, NSC 
PO  = POV(  I ) 

RF  = RFV(  I I 

XS=XV«PC*AV 

ZS=ZV+RC*CV 

XF=XV+RF*AV 

ZF=ZV«PF*CY 

IXS= ( XS+  1 . )*CX 

IZS=(ZS-*1.)*CX 

IXF  = (XF*1.  )*CX 

IZF  = (ZF  + 1.  )*CX 

CALL  GPUF ( IK, ICC,  IXS ,C ) 

CALL  GPLr(IK«l,llC,IZS,0) 

ICX=IXF-IXS 
ICZ=  IZF-  IZS 

CALL  Gt>U  I ( IK42,  S3,  ICX,  IDZ  ) 

IK  = IK4  T 

20  CONTINUE 
RETURN 
FNC 

/FORT 

C 

c 

C CFECKS  ECGE  AGAINST  VIEWING  BCX  SITE 
SUEPCLTINE  CEKS(AN,PN,CN,RC,RF  ) 
CCRPCN/LINE/AV,PV,CV,XV,YY,ZV 
XN=-AN 


r»  o o >•  o o n 


¥^*-eN 

^N=-C^ 

Ca*V*AN«PV*eN-*CV»CN 
rF(C.£G.C. » GC  TO  5 

BI=< (XN-XV)AflN+(YN-YV)*RN+(ZK-7V)*CN)/Q 
IF (C  .L€  .C  . ) CC  TO  3 

C FCGE  CIRECTEC  INTO  VlEkRCX  FRC><  CLTSIDE 
IFIPl.GT.RC)  RC=KI 
PETUPN 
3 CCMtNUF 

C ECGE  CIRECTEC  CLT  CF  VIEtaROX  FBCR  IN51CF 
IF(Rl.LT.RF)  RF=KI 
PETLRN 
5 CCNTINUE 

C FCGE  PARALLEl  TC  SICE 

IF(ARS(XV).GT.l..CR.AflS(ZVl.GT.l.)  RF=RO-i. 

PETLRN 

ENC 

/FORT 


ORCER  ECGE  IMERSECTICN  POINTS  ACCCROING  TO  LOR  AND  FIGE  VALUES 
SLPOCLTINE  CRCL\(RltR?,RC,RF) 

RC  = R 1 


RF  = R2 

IF(R1.LT.R2) 

RC  = P2 

RF='P  1 

RETURN 

ENC 


FORT 


RETURN 


CFECKEC  PARTIALLY  BLUCKEC  LINE  FOR  YISIRLE  S>-GPENTS 
SLERCLTINE  CKLIN(RI.R2I 
CCPPCN/CRAV./NSC,RY,RCV(1C),RFV(10) 

CIPFNSICN  RCdO.RFIiC) 

IF ( R 1 .GT  .C  . .AND  .R2  .LT.RV  I GC  TC  20 
C VERTEX  CE  TEST  ECGE  bEHiND  CONSTRAINT  SURFACE 
IF1P2.GT.BV)  GC  TC  5 
C LCW  ENC  CF  FCGE  FICHEN 
RVC=«2 
RVE=RV 
CC  TC  1C 

C HGF  ENC  CF  FCGE  HICOEN 
5 CONTINUE 
RVC  = C . 

RVF=R1 
1C  CONTINUE 

CC  12  dl.NSC 
RC(  I l=RCV(  I ) 

RF(  I )=RFv(I  ) 

IF  (RFV(  I ) .LE  .RVC  ) RF(l)=-.l. 

IF(ocvm.LT.RVC)  RC(I)=RVC 
IF (RFV( I)  .GT  .RVF  ) RF(1)=RVF 
IF (ROV(  I ) .GT  .RVF  ) RF ( I)=-l. 

12  CCNT  INUE 
GC  TC  3C 

C CONSTRAINT  SURFACE  SEPARATES  ECGE  INTO  TWO  VISIBLE  ENDS 
20  CONTINUF 


46 


BVC  = H I 
RVF  = R2 
K*C 

CC  22  1=1, NSC 

KsK  4 1 

RC(K  } = I ) 

HF (K  )=RFV(  I 1 

IF(RFV([».CT.RVC.«NC.RfV(II.LT.RVF)  Rf(K|=RVC 
IF(RCV(II.Gl.RVC.aNr.RFV(I).L1.aVF)  HF(K)=-1, 

IF (RFV( 1) .n .RVF . 4NC.R0V ( I ) . LI .RVO)  GO  TO  27 

IFIRFVIIt.GI.RVF  .ANG  .ROV  (D.GT.RVO.ANn.ROvm.LT.RVFI  RO(Kl=RVF 
22  CCNT INUF 
NSC=K 
GO  TC  3C 
27  CONT  INC'C 
RF (K  )=PVC 

K=K4l 

RO(K  )=RVF 
RF(K  l=RF\M  I I 
GC  TC  22 

C ARRANGE  RERAIMNG  VISIBLE  SEGKENTS 
30  CCNTIMjF 
NK=C 

CO  32  1 = 1, NSC 

IFIRFI  I I.LE.C.)  GC  TC  32 

NK=NK4l 

POVCNK  )=RCI  I ) 

RFVI NK  ) = RF  t n 
32  CCNTINLF 
NSC=NK 
RETLRN 
ENC 


C FUAT  PLATE  CANCPV  CATA 


c 

NC. 

VERTICES 

, ANC 

SURFACES- 

-BCCV 

f ! 

constraint. 

ANC 

TRANSPARENT 

SC 

30 

23 

■a 

c 

NC. 

VERT ICFS 

P€R 

SURFACE 

A 

3 

3 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

c 

5 

5 

5 

c 

StRFACF 

VERT  ICE 

j 

IN  rotational 

CPCER- 

ALU 

1ST 

VERTICES, 

ALL 

2NC 

, ETC 

• 

1 

1 

2 

ic 

7 

13 

lA 

17 

IF 

22 

22 

2A 

23 

25 

29 

3C 

33 

3A 

5 

37 

39 

A1 

A3 

AS 

A7 

AS 

3E 

39 

AO 

39 

2 

5 

A 

12 

8 

lA 

13 

18 

17 

2A 

21 

23 

21 

26 

30 

29 

3A 

33 

t 

30 

AC 

A2 

AA 

5C 

AP 

At 

t 

A7 

AP 

A1 

j 

3 

e 

11 

10 

1 c 

16 

IS 

2C 

25 

23 

26 

28 

28 

31 

32 

35 

36 

je 

AO 

A2 

AA 

A6 

t 

5C 

A8 

5C 

A9 

A6 

A3 

5 

C 

C 

9 

9 

It 

15 

20 

IS 

27 

2A 

25 

26 

27 

32 

31 

36 

35 

37 

39 

A1 

A3 

A5 

C 

A9 

A7 

AP 

5 

AA 

A5 

C 

0 

C 

C 

C 

C 

C 

C 

C 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C 

0 

c 

C 

C 

C 

C 

C 

AC 

37 

A2 

A7 

C X-PCSITICN  ALL  VERTICES 


7.23 

-7.23 

1A.A7 

13.C79 

-13.079 

5. 

-5. 

5. 

-5. 

5. 

16. 

16. 

16. 

16. 

-16. 

-16. 

-16. 

-16. 

-lA. 

lA. 

-5.5 

5.5 

5.5 

-5.5 

lA. 

-lA. 

16. 

16. 

16. 

16. 

-16. 

-16. 

-16. 

-16. 

1A.A7 

-1A.A7 

16.9AA 

-16.9AA 

16.SA A 

-16.SAA 

16.039 

-16.C39 

1A.366 

-1A.366 

13.273 

-13.273 

12. cot 

-12.CC6 

. 

• 

. 

• 

. 

• 

C Y-PCSITICN  ALL  VERTICES 


25. C6 

25. C6 

61.25 

61.25 

56.769 

56.769 

67.769 

67.769 

67.76S 

67.769 

67.769 

67.769 

85.928 

85.928 

99.928 

99.928 

85.925 

85.925 

99.928 

99.926 

99.928 

99.928 

99.928 

99.928 

113. S28 

113.926 

113.928 

113.926 

133.7 

133.7 

137.7 

137.7 

133.7 

133.7 

137.7 

137.7 

59.06 

59.06 

99.928 

99.928 

154. 6C2 

154. 602 

153.609 

153.609 

152.257 

152.257 

108.032 

108.032 

71.5 

71.5 

• 

• 

• 

• 

• 

. 

Z-PCSITION  ALL  VtRTICES 

63.629 

63.829 

60. 

60. 

76.018 

76.018 

76.018 

76.018 

66.C18 

86.C16 

76.018 

76.016 

71.717 

77.717 

77.717 

71.717 

71.717 

77.717 

77.717 

71.717 

75.717 

75.717 

90.467 

90.467 

92.467 

92.467 

75.717 

75.717 

80.3 

84.6 

84.6 

80.3 

60.3 

84.6 

64.6 

80.3 

67.829 

67.829 

75.717 

75.717 

86.271 

86.271 

93.726 

93.726 

105.96 

105.96 

105.96 

105.96 

92.4 

92  .4 

. 

• 

• 

• 

• 

• 

P ILC7-EYE 

PCS  17  ICN 

0. 

137.7 

95.67 

• 

• 

• 

• 

• 

